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Clarification  regarding  this  year  progress 


My  efforts  were  focused  on  two  main  sub-aims  during  that  year:  1.  studying  the  transcriptional  events 
that  occur  in  cells  upon  ATRX  suppression,  and  analysis  of  the  ALT+  gene  signature  that  I  generated.  This 
effort  is  designed  to  achieve  the  goals  of  Task  2  as  it  appears  in  the  SOW.  2.  Searching  for  positive 
controls  for  suppression  of  ALT  activity.  This  effort  is  relevant  for  Tasks  1  and  3. 


Studying  ALT-specific  transcriptome 

I  analyzed  the  two  gene  expression  datasets  described  in  the  report  to  generate  two  gene  signatures.  I 
used  these  signatures  to  identify  pathways  and  modules  that  are  involved  in  ALT  activity  in  an  attempt 
to  identify  potential  targets  in  ALT-positive  cancers  and  in  order  to  achieve  better  understanding  of  the 
cellular  status  associated  with  ALT  activity.  A  major  finding  of  these  analyses  is  that  HER2  inhibition  as 
well  as  CDK2,  and  HDAC  inhibition  can  reverse  the  ALT+  gene  signature.  This  finding  has  a  therapeutic 
potential  and  I  am  in  the  process  of  validating  it  and  characterizing  it  in  order  to  confirm  its  relevance. 
Similarly,  a  gene-set  enrichment  analysis  identified  an  enrichment  of  DNA  repair  related  genes  in  ALT+ 
cells.  Since  ALT  activity  is  likely  to  increase  the  burden  of  DNA  repair  requirement  in  the  cells,  it  is 
possible  that  the  increased  expression  of  DNA  repair  genes  is  essential  for  the  survival  of  ALT+  cells,  and 
is  therefore  a  potential  target.  The  outcome  of  suppressing  DNA  repair  mechanism  in  ALT+  cells  is 
among  the  questions  I  am  studying  now. 

This  research  was  included  in  the  report. 


Searching  for  positive  controls  for  suppression  of  ALT  activity 

In  parallel  to  the  analysis  of  the  gene  expression  data,  I  am  also  putting  an  effort  toward  identifying  new 
positive  controls  for  Task  3.  As  was  included  in  the  SOW  and  in  previous  reports,  I  tried  to  use  genes  that 
are  known  to  be  involved  in  ALT  activity  as  positive  controls  for  Task  3.  I  obtained  multiple  shRNAs  for 
each  of  the  genes,  however,  as  shown  in  Figure  1  of  the  report,  I  could  not  achieve  a  dynamic  range  that 
is  sufficient  for  an  HTS  screen.  Therefore,  I  am  attempting  to  identify  better  positive  controls  for  ALT 
activity  suppression.  Part  of  this  search  uses  the  gene  expression  data  that  I  generated  in  Task  2  to  guide 
the  selection  of  new  potential  positive  controls  by  testing  genes  whose  expression  in  ALT+  cells  is 
significantly  higher  than  in  non-ALT  cells.  Identifying  robust  positive  controls  with  a  wide  dynamic  range 


is  essential  for  the  success  of  Task  3,  since  it  will  determine  the  threshold  of  calling  genes  as  hits  in  the 
screen,  and  is  therefore  will  determine  the  specificity  (false  positive  rate)  and  sensitivity  (false  negative 
rate)  of  the  screen.  Importantly,  the  genes  that  were  planned  to  be  used  as  positive  controls  in  Task  3 
are  also  the  genes  that  I  planned  to  use  to  suppress  ALT  activity  for  Task  1.  According  to  the  C-circle 
assay,  the  suppression  of  these  genes  resulted  in  only  partial  suppression  of  ALT  activity.  In  order  to 
reliably  identify  proteins  associated  with  PML  and  SplOO  in  ALT+  cells  by  mass-spectrometry  analysis,  it 
is  essential  to  achieve  a  strong  suppression  of  ALT  activity  in  the  perturbed  samples.  As  described  above, 
I  am  using  the  gene  expression  data  to  search  for  genes  whose  perturbation  leads  to  stronger 
suppression  of  ALT  activity.  I  expect  to  proceed  with  both  Task  1  and  Task  3  as  planned  as  soon  as  I 
resolve  this  issue. 

In  addition,  as  I  have  shown  in  previous  reports,  I  adjusted  the  C-circle  assay  to  allow  quantitative  and 
more  sensitive  detection  of  C-circle  abundance.  However,  the  throughput  of  the  assay  was  not  optimal 
for  a  screening  experiment.  I  am  working  on  further  improving  the  assay  and  adjusting  it  for  higher 
throughput  experiments.  The  main  strategy  is  to  enable  detection  of  C-circles  by  qPCR  rather  than  by 
probe  hybridization.  This  optimization  of  the  assay  used  as  a  readout  in  a  HTS  screen  as  planned  in  Task 
3  is  important  for  the  success  of  the  screen,  and  therefore  I  am  investing  in  it. 


This  clarification  regarding  Task  1  and  3  was  added  to  the  revised  report. 
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Introduction 

The  capacity  to  maintain  elongated  telomeres  is  an  essential  feature  of  tumor  cells  enabling 
their  unlimited  replicative  potential.  Approximately  85-90%  of  human  cancers  maintain  telomeres 
through  activation  of  telomerase.  However,  several  cancer  types  harness  telomerase-independent 
mechanisms  referred  to  as  Alternative  Lengthening  of  Telomeres  (ALT).  Importantly,  inhibition 
of  telomerase  in  immortalized  premalignant  cells  and  in  cancer  cells  induces  a  switch  from 
telomerase  activation  to  ALT,  enabling  the  cells  to  escape  therapy  (i,  2). 

The  phenotypic  characteristics  of  ALT-i-  cells  include  long  and  heterogeneous  telomeric 
DNA,  ALT- associated  PML  Bodies  (APBs),  high  level  of  recombination  at  telomeres,  and 
abundant  extrachromosomal  telomeric  DNA,  especially  single  strand  C-rich  circles  (C-circles)(3). 
Currently,  there  is  limited  molecular  understanding  of  the  mechanisms  that  induce  ALT  and  that 
are  involved  in  ALT  activity. 

Importantly,  ALT  activity  has  been  reported  in  breast  carcinomas,  where  it  was  correlated 
with  HER2  positivity  and  with  poor  prognosis  (4).  Furthermore,  a  recent  report  showed  that  ALT 
activity  and  hTERT  expression  can  co-exist  in  the  same  cells  in  patient  breast  tumor  smaples  (5). 
Hence,  understanding  the  mechanisms  that  regulate  and  facilitate  AET  activity  can  promote  the 
development  of  novel  treatments  for  telomerase-independent  breast  cancer.  Such  treatment  has  an 
attractive  potential  in  helping  patients  with  AET-r-poor  prognosis  breast  cancer  and  as  a 
combination  therapy  for  patients  treated  with  telomerase  inhibitors.  It  is  expected  that  inhibition 
of  AET  activity  will  have  minimal  non-desired  effects  on  the  normal  cell  population  that  typically 
has  no  AET  activity. 

The  aim  of  this  project  is  to  characterize  the  mechanisms  that  mediate  and  regulate  AET 
activity  in  breast  cancer,  by  utilizing  proteomic  and  genomic  approaches  to  describe  and  perturb 
AET-i-  cells. 
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Overall  project  summary 

Calibrating  C-circle  assay  to  allow  robust  quantification 

In  order  to  allow  accurate  quantification  of  C-circles  abundance,  the  detection  method  was 
modified  to  enable  digital  acquisition  of  the  signal  with  high  sensitivity.  Using  LiCor  Odessey 
imaging  system  along  with  Image  Studio  software  allows  automatic  quantification  of  the  dot-blot 
labeled  with  IRdye-conjugated  probes.  This  adjustment  improved  the  sensitivity  of  the  assay 
enabling  detection  of  C-circles  in  only  20ng  of  genomic  DNA.  In  addition,  converting  to  a 
fluorescence-labeling  method  allows  safer  processing  of  multiple  samples  at  the  same  time. 

An  Alu  repeat  probe  was  used  in  parallel  to  the  C-circle  probe  to  normalize  the  C-circle  signal 
to  the  input  amount  of  genomic  DNA,  thus  enabling  better  quantitation  of  ALT  activity. 

Testing  these  adjustments  on  ALT-i-  cells,  known  to  have  variable  level  of  C-circles, 
demonstrated  a  wide  range  of  sensitivity  and  very  low  background  from  non- ALT  cells  (Figure 
lA).  Additionally,  I  tested  4  shRNAs  that  target  proteins  known  to  be  involved  in  ALT  activity 
{6-8).  In  contrast  to  LacZ-specific  shRNA  that  had  no  effect  on  C-circle  abundance,  targeting 
RAD51,  MUS81,  SMC5,  and  RAD50  reduced  the  C-circle  signal,  suggesting  a  reduction  in  ALT- 
related  recombination  activity  (Figure  IB). 

Therefore,  the  adjustments  I  introduced  to  the  C-circle  assay  retain  the  sensitivity  of  the  assay 
and  importantly,  improve  quantification  and  assay  handling. 

Study  ATRX  and  DAXX  involvement  in  ALT 

A  study  by  Heaphy  et  al.  (9)  found  that  pancreatic  neuroendocrine  tumors  that  exhibit  features 
of  ALT  activity,  also  have  alternations  in  ATRX  and  DAXX.  It  was  found  that  these  tumors  had 
either  ATRX  or  DAXX  mutations  or  lacked  their  nuclear  expression  at  the  protein  level.  A  similar 
correlation  between  ATRX  and  ALT  activity  was  also  found  in  pediatric  glioblastomas  {10), 
liposarcomas  {11).  A  recent  report  also  showed  that  overexpression  of  ATRX  in  ALT-i-  fibroblasts 
can  suppressed  ALT  activity  {12).  DAXX  and  ATRX  are  histone  chaperones  associated  with 
replication-independent  deposition  of  Histone  H3.3  {13,  14).  ATRX  was  shown  to  be  localized  at 
telomeres  and  to  associate  with  structural  maintenance  of  telomeres  integrity  {15).  To  test  potential 
involvement  of  ATRX  in  breast  cancer,  we  analyzed  its  status  across  the  TCGA  data-set  for  breast 
invasive  carcinoma,  and  found  that  it  is  mutated  in  2%  of  the  tumors  sequenced  by  TCGA  (cBio 
portal  analysis). 

To  further  study  ATRX  involvement  in  ALT  activation,  we  first  tested  its  expression  in  a  panel 
of  22  ALT-I-  cell-lines,  showing  that  in  16  cell-lines  ATRX  expression  was  undetectable  (Figure 
lA  in  {16)).  DAXX  expression  was  reduced  in  only  4  of  the  tested  ALT-i-  cell-lines.  In  3  cell-lines 
that  expressed  ATRX  protein,  it  did  not  exhibit  a  punctuated  nuclear  localization,  suggesting  loss 
of  functional  ATRX  (Figure  IB  in  {16)). 

The  observed  association  between  ALT  activity  and  loss  of  functional  ATRX  encouraged  me 
to  test  whether  knock-down  of  ATRX  or  DAXX  in  non-immortalized  cells  will  promote 
immortalization  through  activation  of  ALT.  To  this  end,  I  used  BJ  cells  stably  expressing  SV40- 
Early  Region  non-immortalized  cell-line.  Control  cells,  infected  with  shRNA  targeting  LacZ,  did 
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not  undergo  immortalization,  but  rather  exhibited  a  gradual  decrease  in  population  doubling  rate, 
and  eventually  stopped  proliferating  at  accumulated  population  doubling  66  (day  73),  typical  to 
these  cells  (Figure  2A). 

To  test  the  hypothesis  that  loss  of  ATRX  or  DAXX  are  involved  in  ALT-mediated 
immortalization,  I  introduced  shRNAs  to  knock-down  ATRX  and  DAXX  in  BJ  cells,  and  the  cells 
were  passaged  in  culture  for  110  days.  Knock-down  of  either  ATRX  or  DAXX  did  not  facilitate  a 
rapid  immortalization  typical  to  hTERT  overexpression.  However,  after  58  and  65  days  in  culture, 
two  immortal  clones  emerged  in  shATRX589-  and  shATRX590-infected  plates,  respectively. 
While  the  population  doubling  rate  before  the  emergence  of  the  two  clones  was  similar  to  that  of 
control  (shLacZ1650-infected)  cells,  this  rate  significantly  increased  to  become  similar  to  the  rate 
of  hTERT-immortalized  cells,  indicating  a  complete  recovery  of  population  crisis  (Eigure  2B, 
compare  week  7  to  week  15). 

Investigating  the  mechanism  by  which  the  emerged  two  colonies  underwent  immortalization, 
I  found  that  the  expression  of  hTERT  significantly  increased  in  the  two  immortalized  shATRX 
colonies  (Eigure  3A).  This  result  suggested  that  the  mechanism  of  telomere  maintenance  in  the 
two  shATRX  immortalized  colonies  is  telomerase-mediated.  To  verify  that  AET  activity  is  not 
involved,  I  used  the  C-circle  assay.  In  accordance  with  telomerase-mediated  telomere 
maintenance,  there  was  no  evidence  for  activation  of  AET  in  the  two  shATRX  immortalized 
colonies  (Eigure  3B). 

Another  feature  related  to  telomere  maintenance,  is  the  expression  of  TERRA  (telomeric 
repeat-containing  RNA),  that  are  transcribed  from  subtelomeric  and  telomeric  regions.  TERRA 
are  associated  with  telomeres  forming  a  heterochromatin  complex,  which  was  suggested  to 
negatively  regulate  telomerase  activity  in  a  telomere-length-dependent  manner  (77).  In  addition, 
in  tumor  cells,  telomerase  activity  correlates  with  lower  TERRA  expression,  compared  to  AET- 
mediated  telomere  maintenance  (18).  Testing  the  expression  of  several  TERRA  transcribed  from 
different  chromosomes,  showed  that  at  early  population  doubling  TERRA  expression  in  control 
cells  and  shATRX-infected  cells  is  similar  (Eigure  3C,  left  panel).  However,  upon  immortalization, 
the  two  shATRX  colonies  exhibited  a  significantly  reduced  expression  of  TERRA  (Eigure  3C, 
right  panel),  in  accordance  with  activation  of  functional  telomerase. 

To  summarize,  knock-down  of  either  ATRX  or  DAXX  did  not  facilitate  rapid  immortalization 
and  AET  activation  in  cells.  Eack  of  ATRX  activity  was  shown  before  to  be  associated  with  modest 
level  of  telomere  dysfunction  (75).  A  recent  report  by  Clynes  et.al.  also  showed  that  ATRX 
suppression  is  not  sufficient  to  induce  AET  activity  (79).  Moreover,  this  report  demonstrated  that 
ATRX  suppression  leads  to  DNA  damage  response  in  telomeres.  To  survive,  the  cells  are  forced 
to  activate  a  mechanism  that  elongate  and  maintain  telomeres.  This  resolution  is  typically  a  rare 
event  which  occurs  at  a  single  cell  allowing  a  significant  increase  in  proliferation  rate  that 
eventually  result  in  the  emergence  of  an  immortalized  colony  that  will  rapidly  take-over,  the 
otherwise  arrested,  population.  Though  the  cells  I  used  resolved  this  by  activation  of  telomerase, 
it  is  nevertheless  possible  that  under  other  cellular  contexts,  this  could  be  resolved  by  activation 
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of  ALT,  as  is  indicated  by  the  correlation  of  ALT  occurrence  with  ATRX  mutations  and  lack  of 
expression. 


Unbiased  analysis  of  transcriptional  profiles  of  ALT+  and  non-ALT  cells 

In  order  to  characterize  the  unique  transcriptional  profile  of  ALT+  cells,  I  measured  the  steady- 
state  global  transcriptional  program  of  16  ALT-i-  and  9  non-ALT  cell-lines  using  the  Affymetrix 
microarray  HTU133A  platform. 

In  order  to  visualize  the  variance  in  the  data,  I  used  Principal  Component  Analysis  (PCA), 
which  allows  the  unbiased  identification  of  major  sub-groups  of  cell-lines  within  the  complete  data 
set  that  exhibit  higher  correlation  to  each  other.  As  shown  in  Figure  4A,  the  first  two  PCAs  with 
the  highest  variance  clustered  the  cell-lines  into  two  major  groups  according  to  their  ALT  activity. 
Importantly,  this  analysis  shows  that  the  ALT-i-  cell-lines  correlate  better  with  each  other  than  with 
non-ALT  cell-lines.  Furthermore,  this  analysis  also  indicates  that  the  most  significant  feature  that 
distinguishes  the  samples  is  their  ALT  activity,  since  the  two  most  variant  PCAs  were  able  to 
distinguish  the  ALT-i-  from  the  non-ALT  cell-lines. 

To  explore  the  unique  transcriptional  profile  related  to  ALT  activity,  I  first  identified  the  genes 
that  exhibit  variability  within  the  data  using  GenePattern  (20).  A  total  of  9,909  probes  exhibited 
over  two-fold  change  and  a  minimal  delta  of  50  across  cell-lines  (pre-normalized  data), 
independent  of  ALT-activity  status.  In  accordance  with  the  PCA  analysis  (Figure  4A),  a  correlation 
matrix  of  all  the  samples  showed  two  major  clusters  distinguishing  cell-lines  by  their  ALT  activity 
(Figure  4B).  The  correlation  coefficients  presented  here  are  affected  by  cell-line  specific 
expression  patterns  as  well  as  by  ALT-related  patterns  since  all  the  genes  that  exhibited  variability 
across  all  samples  were  included  in  computing  this  correlation  matrix.  This  suggests  that  the 
contribution  of  ALT-related  expression  patterns  is  significant  enough  to  lead  an  unbiased 
clustering  of  the  measured  cell  lines. 

Studying  the  differential  expression  pattern  of  ALT-i-  cells 

To  facilitate  better  understanding  of  the  molecular  mechanisms  that  are  related  to  ALT 
activity,  I  used  the  GenePattern  Comparative  Marker  Selection  module  to  identify  genes  that 
exhibit  differential  expression  between  ALT-i-  and  non-ALT  cell-lines  (2i).  The  genes  were  then 
ranked-ordered  by  this  score  (Figure  5A).  Based  on  this  analysis,  I  identified  5022  genes  whose 
differential  expression  significantly  (q-value<0.01)  correlates  with  ALT  activity,  of  which  2147 
genes  had  a  q-value<0.001  (corresponding  to  FDR<0.1%;  Figure  5B). 

Investigating  the  list  of  differentially  expressed  genes  revealed  that  several  known  ALT- 
related  genes  were  indeed  up-regulated  in  ALT-i-  cells.  Members  of  the  MRN  complex,  which  is 
essential  for  telomere  elongation  in  ALT-i-  cell-lines  (6),  MREll,  RAD50,  and  XRCCl,  were 
significantly  up-regulated  in  ALT-i-  compared  to  non-ALT  cell-lines  (Figure  5C-E).  Additional 
genes  associated  with  AET  activity,  such  as  SMC5/6,  EENl,  ASEIB,  RECQE,  MUS81,  TREEl/2, 
EANCE,  and  SPlOO,  were  also  significantly  up-regulated  in  the  AET-i-  cell-lines. 
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To  further  study  the  unique  gene  expression  pattern  of  ALT+  cell-lines,  I  took  an  unbiased 
approach  to  analyze  the  gene  signature  identified  above  by  using  GSEA  (Gene  Set  Enrichment 
Analysis)  tool.  Eirst,  I  analyzed  the  enrichment  of  the  AET-related  signature  genes  against  all  Gene 
Ontology  (GO)  gene  sets.  Each  of  the  sets  in  this  collection  consists  of  genes  that  are  annotated  by 
the  same  GO  term.  24  out  of  the  1025  gene  sets  included  in  this  collection  were  significantly 
enriched  for  in  the  AET-i-  cell-lines  with  a  nominal  p-value  lower  than  1%.  Several  of  the  enriched 
gene  sets  were  associated  with  DNA  replication,  mitosis,  chromosome  segregation,  and  spindle 
regulation  (Table  1).  Interestingly,  one  of  the  top  enriched  gene  sets  identified  here  is  termed  ‘ATP 
DEPENDENT  DNA  HELICASE  ACTIVITY’  (Eigure  6A),  which  includes  genes  whose  activity 
relates  to  the  MRN  complex  as  well  as  to  BEM  and  WRN. 

A  second  gene  set  collection  I  used  consists  of  curated  gene  sets  that  are  based  on  data 
from  the  literature,  pathway  databases  (such  as  KEGG  and  REACTOM),  and  perturbation  datasets. 
153  out  of  the  3767  gene  sets  included  in  this  collection  were  significantly  enriched  for  in  the 
AET-i-  cell-lines  with  a  nominal  p-value  lower  than  1%  (Table  2).  Several  of  the  enriched  gene  sets 
identified  in  this  analysis  are  related  to  DNA  synthesis  and  more  specifically,  to  the  synthesis  of 
the  lagging  strand  (Eigure  6B).  Confirming  the  relevance  of  this  gene  set  to  AET  activity,  the  set 
includes  EENl,  which  is  essential  for  stabilizing  telomeres  in  AET-i-  cells  but  not  in  normal  or 
telomerase-active  cells  (22).  Additional  enriched  gene  sets  that  could  potentially  be  relevant  to 
AET-t  cell-lines  are  sets  related  to  DNA  repair  and  maintenance  of  stability  (Eigure  6C).  AET 
activity  involves  extensive  recombination  at  the  chromosomal  ends,  and  thus  could  increase  the 
overall  burden  of  DNA  repair  in  the  cell.  Therefore,  over-expression  of  repair  genes  could  be 
essential  for  maintaining  DNA  stability  in  AET-i-  cells.  Of  note,  several  of  the  genes  included  in 
the  set,  such  as  AURKA,  are  already  being  tested  as  targets  for  targeted  therapy  in  breast  cancer, 
suggesting  AET-i-  breast  tumors  could  benefit  from  these  drugs  (23,  24). 

Einally,  a  gene  set  consisting  of  genes  encoding  kinases  that  are  differentially  expressed  in 
basal  breast  cancer  was  enriched  for  in  the  AET-i-  cell-lines  (Eigure  6D)  (25).  Most  of  these  kinases 
are  associated  with  cell-cycle  regulation,  such  as  AURKA,  AURKB,  BUBl,  and  CDKl.  PEKl 
was  reported  to  contribute  to  AET  activity  by  phosphorylating  TREFl,  thus  promoting  it’s  binding 
and  stabilization  of  telomeres  (26).  Eikewise,  CDKl  regulates  the  replication  of  telomeric  DNA 
in  telomerase-negative  cells  (27),  and  also  phosphorylates  TREEl  to  promote  its  recruitment  to 
sites  of  DNA  damage  (28). 

The  differential  expression  of  these  genes  between  AET-i-  and  non-AET  cell-lines  could 
potentially  be  utilized  to  gain  therapeutic  specificity,  however,  further  investigation  of  the 
contribution  of  these  kinases  to  AET  activity  is  necessary. 

To  summarize,  AET-i-  cell-lines  exhibit  a  differential  expression  pattern  as  compared  to  non- 
AET  cell-lines.  This  AET-related  expression  signature  includes  several  known  AET-related  genes, 
which  signifies  its  relevance  to  understanding  the  mechanisms  that  mediate  and  regulate  AET 
activity.  Using  statistical  analysis  tools,  I  identified  several  sets  of  genes  that  are  differentially 
expressed  in  AET-i-  cell-lines  and  are  potentially  involved  in  AET  activity  and  the  survival  of 
AET-t  cells. 
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Predicting  sensitivity  to  drugs  using  the  ALT+  gene  signature 

The  Connectivity  Map  dataset  contains  over  6000  gene  expression  profiles  from  human  cells 
treated  with  small  molecule  compounds  (29).  The  connectivity  Map  tool  compares  a  query  gene 
signature  of  interest,  consisting  of  up-  and  down-regulated  genes,  to  all  the  gene  signatures  in  its 
dataset.  The  connection  between  the  query  signature  and  each  of  the  dataset  signatures  is  calculated 
using  rank-based  pattern  matching  and  scored  according  to  Kolmogorov-Smirnov  statistics. 

Using  the  ALT-i-  gene  expression  signature  generated  above  as  a  query  signature,  I  generated 
a  list  of  small  molecule  compounds  whose  gene  signatures  are  negatively  connected  to  the  ALT-i- 
signature  (Table  3).  Potentially,  negatively  connected  compounds  reverse  the  given  query 
signature,  suggesting  that  few  of  the  identified  compounds  could  convert  the  ALT-i-  state  of  gene 
expression  into  non- ALT  cellular  state.  Interestingly,  3  of  the  top  20  compounds  that  had  a  strong 
negative  connectivity  to  the  ALT-i-  signature  have  known  relevance  to  breast  cancer;  The  strongest 
negative  connection  identified  is  the  perturbation  of  cells  with  the  HER2  inhibitor  AG-825.  This 
connectivity  suggests  that  treating  HER2  positive,  AET-dependent  breast  tumors  with  an  HER2 
inhibitor  would  affect  the  tumor  cells  both  by  inhibiting  the  oncogenic  HER2  activity,  and  by 
reverting  the  AET-specific  signature.  Eikewise,  the  HD  AC  inhibitor  MS-275  (Entinostat)  and  the 
CDK2  inhibitor  GW-8510  are  also  negatively  connected  to  the  AET-i-  gene  signature.  Both 
inhibitors  have  been  tested  in  breast  cancer  (30,  31),  and  potentially  could  have  greater  efficiency 
in  AET-i-  breast  tumors  due  to  their  ability  to  revert  the  AET-i-  cellular  state. 


Analysis  of  gene  expression  patterns  upon  ATRX  knock-down 

Several  reports  indicate  the  involvement  of  ATRX  in  AET  activity  and  in  telomere 
maintenance.  Since  ATRX  is  a  chromatin  remodeler  it  is  expected  that  suppression  of  ATRX  will 
affect  the  expression  of  multiple  genes,  including  genes  that  regulate  AET  activity.  Although 
knocking-down  ATRX  in  BJ  cells  resulted  in  immortalization  facilitated  by  hTERT  expression 
rather  than  AET  activation  (see  above  and  (16)),  I  expected  that  the  suppression  of  ATRX  would 
still  affect  the  expression  of  genes  that  are  involved  in  AET  activity.  Therefore,  I  profiled  the 
steady-state  transcriptome  of  both  parental  BJ  cells  and  shATRX589-immortalized  cells  by 
RNAseq  using  Illumina  HiSeq  platform.  Reads  were  aligned,  mapped,  and  counted  using  the 
Tuxedo  tools.  The  10,000  highest-expressing  genes  were  considered  for  further  analysis  of 
differential  expression  in  order  to  ensure  that  signals  are  above  noise  (Eigure  7A).  To  identify 
differentially  expressing  genes,  I  calculated  the  fold-change  in  expression  upon  shATRX-mediated 
immortalization.  There  are  745  genes  that  exhibit  more  than  1.5  fold-change  upon  suppression  of 
ATRX,  of  which  233  genes  were  up-regulated  and  512  genes  were  down-regulated  (Eigure  7B). 

To  explore  the  genes  that  are  differentially  expressed  upon  ATRX  suppression,  I  used  GSEA 
to  identify  enrichment  of  gene  sets.  Interestingly,  one  of  the  genes  sets  that  was  significantly 
enriched  for  upon  ATRX  suppression  and  immortalization  is  named  ‘PID_Telomerase_pathway’ 
and  contained  several  AET-related  protein  such  as  TREE  1/2  and  members  of  the  XRCC  complex. 


12 


as  well  as  nibrin  (NBN)  a  member  of  the  MRN  complex  which  is  essential  for  ALT  activity  (6). 
The  up-regulation  of  ALT-associated  genes  upon  ATRX  suppression  confirms  the  potential 
involvement  of  ATRX  in  ALT  activation,  and  importantly  suggests  that  more  ALT-related  genes 
can  be  identified  by  further  mining  the  gene  expression  changes  following  ATRX  suppression. 


Searching  for  positive  controls  for  suppression  of  ALT  activity 

Based  on  the  literature,  I  tried  to  use  genes  that  are  known  to  be  involved  in  ALT  activity  as 
positive  controls  for  determining  ALT  inhibition  in  Aim  3. 1  obtained  multiple  shRNAs  for  each 
of  the  genes,  however,  in  my  hands  and  as  shown  in  Figure  IB,  I  could  not  achieve  a  dynamic 
range  that  is  sufficient  for  an  HTS  screen.  Therefore,  I  am  attempting  to  identify  better  positive 
controls  for  ALT  activity  suppression.  Part  of  this  search  uses  the  gene  expression  data  that  I 
generated  above  to  guide  the  selection  of  new  potential  positive  controls.  I  am  testing  genes 
whose  expression  in  ALT-i-  cells  is  significantly  higher  than  in  non-ALT  cells.  Identifying  robust 
positive  controls  with  a  wide  dynamic  range  is  essential  for  the  success  of  the  screen,  since  it  will 
determine  the  threshold  of  calling  genes  as  hits  in  the  screen,  and  therefore  will  determine  the 
specificity  (false  positive  rate)  and  sensitivity  (false  negative  rate)  of  the  screen. 

Importantly,  characterization  of  the  proteins  associated  with  PML  and  SplOO  (Aim  1)  relays  on 
the  same  positive  control  genes.  According  to  the  C-circle  assay,  the  suppression  of  these  genes 
resulted  in  only  partial  suppression  of  ALT  activity.  In  order  to  reliably  identify  proteins 
associated  with  PML  and  SplOO  in  ALT-i-  cells  by  mass-spectrometry  analysis,  it  is  essential  to 
achieve  a  strong  suppression  of  ALT  activity  in  the  perturbed  samples.  As  described  above,  I  am 
using  the  gene  expression  data  to  search  for  genes  whose  perturbation  can  lead  to  stronger 
suppression  of  ALT  activity. 
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Key  research  accomplishments 

•  Calibrated  a  fluorescence-based  C-circle  assay  to  enable  quantification  and  safe  detection 
of  ALT  activity  with  moderate  throughput. 

•  Tested  the  contribution  of  ATRX  and  DAXX  suppression  to  ALT  activation. 

•  Profiled  gene  expression  patterns  of  16  ALT-i-  cell-lines  as  well  as  9  non- ALT  cell-lines. 

•  Defined  an  ALT-related  gene  expression  signature. 

•  Identified  several  sets  of  genes  that  are  uniquely  expressed  in  ALT-i-  cells  and  are 
potentially  related  to  ALT  activity. 

•  Identified  small  molecule  inhibitors  negatively  connected  to  the  ALT-i-  gene  expression 
profile,  with  potential  efficiency  in  ALT-i-  tumors. 

•  Profiled  changes  in  gene  expression  upon  suppression  of  ATRX. 
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Conclusion 

In  order  to  shed  more  light  on  the  mechanism  that  lead  to  ALT  activation,  I  focused  on  changes 
in  gene  expression  upon  ALT  activation  by  profiling  the  transcriptome  of  16  ALT+  cell-lines  and 
9  non-ALT  cell-lines.  By  applying  unbiased  analysis  strategy,  I  showed  that  ALT  activity  imposes 
a  significant  and  wide  effect  on  global  gene  expression  patterns.  Further  analysis  revealed  an 
ALT-i-  differential  gene  expression  signature  that  contained  multiple  genes  that  are  already  known 
to  be  involved  in  ALT  activity.  Mining  the  ALT-related  expression  signature  using  Gene  Set 
Enrichment  Analysis  identified  several  sets  of  genes  with  potential  relevance  to  ALT  activity. 
Interestingly,  many  of  these  genes  regulate  DNA  replication,  the  cell-cycle,  and  DNA  damage 
repair.  The  up-regulation  of  those  genes  in  ALT-i-  cell-lines  suggests  that  ALT  activity  involves 
re-wiring  of  these  exiting  mechanisms  to  elongate  and  maintain  telomeres  in  the  absence  of 
telomerase.  It  could  also  indicate  potential  dependencies  of  ALT-i-  cells,  such  as  increased  DNA 
damage  repair  activity,  in  line  with  our  previous  observation  of  high  level  of  genomic  instability 
and  frequent  presence  of  DNA  damage  in  ALT-i-  cells  {16).  Querying  the  connectivity  map  dataset 
with  the  ALT-r  signature  identified  few  small  molecule  inhibitors  that  could  have  increased 
efficiency  in  ALT-i-  breast  cancer  tumors. 

ATRX  is  a  chromatin  modulator  that  was  shown  to  be  genetically  and  functionally  altered  in 
ALT-r  tumors.  Never  the  less,  its  suppression  was  not  sufficient  to  induce  immortalization  of  BJ 
cells,  suggesting  that  additional  factors  are  required.  To  study  the  cellular  events  that  occur  upon 
loss  of  ATRX  function,  I  profiled  global  gene  expression  from  BJ  cells  that  were  immortalized 
upon  knock-down  of  ATRX.  GSEA  analysis  identified  a  set  of  genes  that  contains  several  AET- 
related  genes,  thus  enabling  a  focused  search  within  this  gene  set  for  ATRX-depended  AET- 
regulating  genes. 

The  high  level  of  genomic  instability  and  frequent  presence  of  DNA  damage  foci  observed  by 
us  in  AET  cells  {16),  suggests  that  these  cells  will  have  unique  vulnerabilities,  which  are  yet  to  be 
found. 

Together,  or  comprehensive  study  of  the  genomic  features  of  AET  cells  {16)  and  the  broad 
profiling  of  the  AET-related  gene  expression  are  useful  resources  for  further  studying  the 
mechanism  that  regulates  and  facilitates  AET  activity.  Several  of  the  gene  sets  identified  here 
include  known  drug  targets  and  thus  could  have  important  clinical  implications. 
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Publications,  Abstracts  and  Presentations 

Lovejoy  CA,  Li  W,  Reisenweber  S,  Thongthip  S,  Bruno  J,  de  Lange  T,  De  S,  Petrini  JH,  Sung  PA, 
Jasin  M,  Rosenbluh  J,  Zwang  Y,  Weir  BA,  Hatton  C,  Ivanova  E,  Macconaill  L,  Hanna  M,  Hahn 
WC,  Lue  NF,  Reddel  RR,  Jiao  Y,  Kinzler  K,  Vogelstein  B,  Papadopoulos  N,  Meeker  AK, 
Consortium  ALTSC  (2012)  Loss  of  ATRX,  genome  instability,  and  an  altered  DNA  damage 
response  are  hallmarks  of  the  alternative  lengthening  of  telomeres  pathway.  PLoS  genetics  8: 
el 002772 
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Inventions,  Patents  and  Licenses 

Nothing  to  report 
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Reportable  outcomes 

Nothing  to  report 
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Other  achievements 

Nothing  to  report 
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Opportunities  for  training  and  professional  development 

Training: 

Weekly  scientific  updating  meetings  with  mentor  and  group  members  to  discuss  my  progress  as 
well  as  other  projects  conducted  in  the  group.  Since  several  group  members  are  physician- 
scientists  many  of  the  discussions  were  focused  on  the  relevance  of  the  project  and  results  to 
cancer  detection  and  therapy. 

Bi-weekly  ono-on-one  meetings  with  mentor  to  review  results,  design  experimental  strategy,  and 
discuss  progress. 

Professional  development: 

One  day  workshops  on  bioinformatics  tools  given  by  developer  teams  at  the  Broad  Institute: 
GSEA,  Gene-E,  IGV,  Eirehose,  GenePattem. 

A  week-long  workshop  on  gene-expression  analysis  and  relevant  statistical  tools  given  by 
computational-biology  scientists  at  the  Broad  Institute. 

Cancer  Dependencies  working  group  seminars:  monthly  seminars  on  current  research  and  data 
analysis  on  cancer  dependencies,  potential  targets,  and  therapies. 
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A.  GO:  ATP  DEPENDENT  DNA  HELICASE  ACTIVITY 


B.  REACTOME:  PROCESSIVE  SYNTHESIS  ON  THE  LAGGING  STRAND 


C.  KEGG:  MISMATCH  REPAIR 


D.  FINETTI  BREAST  CANCER  KINOME  RED 


Figure  6 
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NAME 

SIZE 

NES 

MITOTIC  SISTER  CHROMATID  SEGREGATION 

11 

-0.77373 

SISTER  CHROMATID  SEGREGATION 

11 

-0.77373 

ATP  DEPENDENT  DNA  HELICASE  ACTIVITY 

8 

-0.74028 

GI  PHASE  OF  MITOTIC  CELL  CYCLE 

8 

-0.71781 

REPLICATION  FORK 

14 

-0.7006 

ORGANELLAR  SMALL  RIBOSOMAL  SUBUNIT 

8 

-0.68236 

MITOCHONDRIAL  SMALL  RIBOSOMAL  SUBUNIT 

8 

-0.68236 

SMALL  RIBOSOMAL  SUBUNIT 

8 

-0.68236 

CHROMOSOME  SEGREGATION 

20 

-0.67529 

REGULATION  OF  MITOSIS 

28 

-0.6676 

SPINDLE  MICROTUBULE 

15 

-0.66266 

MICROTUBULE  MOTOR  ACTIVITY 

12 

-0.65116 

KINETOCHORE 

20 

-0.64197 

DNA  DEPENDENT  ATPASE  ACTIVITY 

17 

-0.63868 

MITOTIC  CELL  CYCLE  CHECKPOINT 

16 

-0.63795 

GI  PHASE 

9 

-0.63567 

DNA  DIRECTED  DNA  POLYMERASE  ACTIVITY 

12 

-0.62905 

MICROTUBULE  ORGANIZING  CENTER  ORGANIZATION  AND  BIOGENESIS 

9 

-0.62795 

DRUG  BINDING 

10 

-0.62713 

DNA  POLYMERASE  ACTIVITY 

14 

-0.62648 

Table  1:  Top  20  GO  gene  sets  enriehed  for  in  ALT+  eell-lines.  Size  indicates  the  number  of  genes 
included  in  each  set.  NES,  Normalized  Enrichment  Score.  Reflecting  the  degree  to  which  a  gene 
set  is  over-represented  at  either  ends  of  the  ranked  list,  normalized  by  permutations  of  the  ranked 
dataset.  A  more  negative  NES  indicates  higher  enrichment  of  the  gene  set  in  the  AET-i- 


differentially  up-regulated  end  of  the  ranked  list. 
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NAME 

SIZE 

NES 

REACTOME  PROCESSIVE  SYNTHESIS  ON  THE  LAGGING  STRAND 

13 

-0.89422 

REACTOME  REMOVAL  OF  THE  FLAP  INTERMEDIATE  FROM  THE  C  STRAND 

9 

-0.89365 

REACTOME  LAGGING  STRAND  SYNTHESIS 

17 

-0.79381 

REACTOME  POL  SWITCHING 

12 

-0.79306 

LY  AGING  MIDDLE  DN 

16 

-0.79292 

KANG  DOXORUBICIN  RESISTANCE  UP 

47 

-0.78178 

AMUNDSON  GAMMA  RADIATION  RESPONSE 

34 

-0.77391 

FINETTI  BREAST  CANCER  KINOME  RED 

14 

-0.76779 

KALMA  E2F1  TARGETS 

11 

-0.75699 

MONTERO  THYROID  CANCER  POOR  SURVIVAL  UP 

10 

-0.75409 

REACTOME  PURINE  RIBONUCLEOSIDE  MONOPHOSPHATE  BIOSYNTHESIS 

9 

-0.75289 

KUMAMOTO  RESPONSE  TO  NUTLIN  3A  DN 

9 

-0.7507 

BURTON  ADIPOGENESIS  PEAK  AT  24HR 

25 

-0.74505 

REACTOME  REPAIR  SYNTHESIS  FOR  GAP  FILLING  BY  DNA  POL  IN  TC  NER 

14 

-0.74457 

KEGG  MISMATCH  REPAIR 

21 

-0.744 

TANG  SENESCENCE  TP53  TARGETS  DN 

45 

-0.7432 

KAMMINGA  EZH2  TARGETS 

34 

-0.73923 

CROONQUIST  NRAS  SIGNALING  DN 

66 

-0.73693 

REACTOME  SLBP  DEPENDENT  PROCESSING  OF  REPLICATION  DEPENDENT 

HISTONE  PRE  MRNAS 

8 

-0.73269 

CROONQUIST  IL6  DEPRIVATION  DN 

80 

-0.7323 

Table  2:  Top  20  curated  gene  sets  enriched  for  in  ALT+  cell-lines. 
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Name 

Enrichment  score 

p-value 

tyrphostin  AG-825 

-0.967 

— 

piperlongumine 

-0.961 

0.0033 

dexverapamil 

-0.959 

— 

MS-275 

-0.946 

0.00634 

celastrol 

-0.946 

— 

menadione 

-0.913 

0.01521 

trazodone 

-0.891 

0.00248 

fisetin 

-0.888 

— 

ebselen 

-0.868 

0.00457 

1 ,4-chrysenequinone 

-0.868 

0.03455 

GW-8510 

-0.867 

0.00062 

cantharidin 

-0.864 

— 

rottlerin 

-0.861 

0.00547 

blebbistatin 

-0.854 

0.04269 

luteolin 

-0.853 

0.00084 

clomipramine 

-0.85 

0.00093 

thioguanosine 

-0.85 

0.00093 

apigenin 

-0.849 

0.00095 

sulconazole 

-0.845 

0.00103 

triflusal 

-0.839 

0.00841 

Table  3:  Top  20  small  compound  perturbations  negatively  connected  to  the  ALT+  gene  expression 
signature.  Enrichment  score  represents  the  strength  of  connectivity. 
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Legend  to  figures 

Figure  1:  Calibration  a  fluorescence-based  C-circle  assay.  A.  Genomic  DNA  was  isolated 
from  the  indicated  cells.  Following  C-circle  assay  reaction,  samples  were  dot-blotted,  and  an 
IRDye-labeled  probe  specific  for  either  C-circle  assay  product  or  Alu  repeats  was  hybridized  to 
the  membranes.  Membranes  were  scanned  using  LiCor  Odessey  imaging  system,  and  dotes 
quantified  using  Image  Studio  software.  Bar-graph  presents  relative  abundance  of  C-circles.  B. 
indicated  shRNAs  were  virally  delivered  into  U20S  cells,  followed  by  2  days  of  puromycin 
selection  and  additional  2  days  propagation  under  normal  culture  conditions.  Genomic  DNA  was 
then  isolated,  and  C-circle  assay  performed  and  quantified  as  in  A. 

Figure  2:  Immortalization  of  BJ  cells  upon  knock-down  of  ATRX  and  DAXX.  A.  BJ  cells 
were  infected  with  LacZ-  ATRX-  or  DAXX-specific  shRNAs  at  population-doubling  (PD)  40. 
Cells  were  selected  for  2  days  with  puromycin,  and  at  least  1^’  cells  for  each  sub-line  were 
maintained  in  culture  for  1 10  days.  Cultures  were  split  at  <80%  confluence,  and  growth  media  was 
re-freshed  every  3  days.  Cells  were  counted  at  each  culture  passaging  and  accumulated  PD  was 
recorded.  B.  Rate  of  population-doubling  per  day  is  plotted  for  the  indicated  sub-lines  at  the 
indicated  times,  corresponding  to  boxed  time-points  in  A. 

Figure  3:  Characterizing  telomere  maintenance  mechanism  of  the  two  immortalized  shATRX 
colonies.  A.  RNA  was  isolated  from  the  indicated  sub-lines  at  weeks  2  and  15  (see  Figure  2A). 
hTERT  expression  was  measured  using  qPCR  with  either  an  internal  primer  set  designed  around 
exons  5  and  6,  or  a  3’UTR  primer  set  designed  around  the  last  exon  and  the  3’UTR  of  the 
endogenous  gene.  Shown  is  fold-change  of  expression  level  at  week  15  compared  to  week  2  for 
the  two  probe  sets.  B.  genomic  DNA  was  isolated  at  week  15  from  the  indicated  sub-lines  (see 
Figure  2A).  Shown  are  dot-blots  of  the  C-circle  probe  and  Alu  probe.  JFCF/T.IF  and  JFCF/T.IJ 
were  used  as  negative  and  positive  controls,  respectively.  C.  RNA  was  isolated  from  the  indicated 
sub-lines  14  days  after  infection  or  after  immortalized  colonies  emerged  in  the  shATRX  sub-lines. 
TERRA-specific  cDNA  was  synthesized  and  TERRA  expressed  from  chromosome  17p  was 
amplified. 

Figure  4:  Analysis  of  gene  expression  profiles  of  AET-i-  and  non-AET  cell-lines  shows 
significant  AET-related  difference  in  global  gene  expression.  A.  PCA  analysis  of  normalized 
expression  data.  The  values  for  the  PCAs  with  the  highest  variance  are  plotted  for  each  cell-line. 
Green,  AET-i-  cell-lines;  Orange,  non-AET  cell-lines.  The  circles  labeled  I  and  II  mark  the  two 
identified  clusters.  B.  Correlation  coefficients  were  calculated  for  all  variable  genes  across  all  cell¬ 
lines  and  plotted  as  a  heatmap.  Red,  positive  correlation;  Blue,  Negative  correlation.  Dendrogram 
above  corresponds  to  the  degree  of  linkage  between  cell-lines. 

Figure  5:  Comparative  Marker  Selection  analysis  defines  an  AET-related  gene  signature.  A. 
Plotted  is  the  expression  of  genes  across  all  cell-lines  ranked  according  to  the  correlation  between 
the  gene’s  expression  and  the  phenotypic  classification  of  the  cell-lines.  Red,  up-regulated;  Blue, 
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down-regulated.  B.  The  expression  of  the  2147  most  significantly  differentially  expressed  genes 
according  to  the  Comparative  Marker  Selection  analysis  is  plotted.  C-E.  Expression  of  known 
ALT-related  genes  is  compared  between  non-ALT  and  ALT-i-  cell-lines.  Mean  and  SD  indicated 
by  bar,  p-value  for  differential  expression  (T-test)  is  presented. 

Figure  6:  GSEA  reveals  potential  associations  to  AET  activity.  A-D.  Expression  of  genes 
included  in  the  indicated  enriched  gene  sets  across  all  cell-lines.  Heatmap  colors  as  in  Eigure  5A. 

Figure  7:  Gene  expression  profiling  of  cells  immortalized  by  ATRX  knock-down.  A.  Gene 
expression  of  parental  BJ  cells  and  BJ-shATRX589  immortalized  cells.  Genes  ranked-ordered 
according  to  their  expression  fold-change  following  suppression  of  ATRX.  Red,  up-regulated; 
Blue,  down-regulated.  B.  Expression  of  genes  exhibiting  over  1.5  fold-change  upon  ATRX  knock¬ 
down. 
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Abstract 

The  Alternative  Lengthening  of  Telomeres  (ALT)  pathway  is  a  telomerase-independent  pathway  for  telomere  maintenance 
that  is  active  in  a  significant  subset  of  human  cancers  and  in  vitro  immortalized  cell  lines.  ALT  is  thought  to  involve 
templated  extension  of  telomeres  through  homologous  recombination,  but  the  genetic  or  epigenetic  changes  that  unleash 
ALT  are  not  known.  Recently,  mutations  in  the  ATRX/DAXX  chromatin  remodeling  complex  and  histone  H3.3  were  found  to 
correlate  with  features  of  ALT  in  pancreatic  neuroendocrine  cancers,  pediatric  glioblastomas,  and  other  tumors  of  the 
central  nervous  system,  suggesting  that  these  mutations  might  contribute  to  the  activation  of  the  ALT  pathway  in  these 
cancers.  We  have  taken  a  comprehensive  approach  to  deciphering  ALT  by  applying  genomic,  molecular  biological,  and  cell 
biological  approaches  to  a  panel  of  22  ALT  cell  lines,  including  cell  lines  derived  in  vitro.  Here  we  show  that  loss  of  ATRX 
protein  and  mutations  in  the  ATRX  gene  are  hallmarks  of  ALT-immortalized  cell  lines.  In  addition,  ALT  is  associated  with 
extensive  genome  rearrangements,  marked  micronucleation,  defects  in  the  G2/M  checkpoint,  and  altered  double-strand 
break  (DSB)  repair.  These  attributes  will  facilitate  the  diagnosis  and  treatment  of  ALT  positive  human  cancers. 
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Introduction 

In  the  absence  of  telomerase  activity,  telomeres  shorten  with  ceU 
division,  ultimately  leading  to  senescence.  Hence,  the  development 
of  human  cancer  is  associated  with  an  active  telomere  mainte¬ 
nance  system  that  provides  an  infinite  source  of  telomeric  DNA  to 
potentiate  immortality.  Although  telomerase  reactivation  is  the 
most  common  mechanism  of  telomeric  repeat  addition  in  cancers, 
a  significant  subset  of  human  tumors  employs  a  telomerase- 
independent  telomere  maintenance  pathway,  referred  to  as  ALT 
[!]■ 

The  emerging  view  is  that  ALT  maintains  telomeres  through 
homology-directed  recombination  (HDR)  [2].  Supporting  this 


view,  ALT  cells  show  an  elevated  frequency  of  sequence  exchanges 
between  telomeres  [2-4],  contain  extrachromosomal  linear  and 
circular  telomeric  DNA  [5-8],  and  often  exhibit  heterogeneously- 
sized  telomeres  [1],  features  consistent  with  hyperactive  HDR. 
The  extrachromosomal  telomeric  DNA  has  been  proposed  to 
serve  as  a  template  for  the  extension  of  telomeres  by  a 
recombination  mechanism  akin  to  Break  Induced  Replication 
(BIR)  [9,10].  In  addition,  ALT  cells  often  carry  altered  PML 
bodies  (ALT-associated  PML  bodies,  or  APBs)  that  contain 
telomeric  DNA  as  well  as  numerous  recombination  factors  [11]. 
Several  proteins  involved  in  recombination  are  known  to  be 
required  for  ALT,  including  the  Mre  1 1  complex,  Mus8 1 ,  and  the 
SMG5/6  sumoylation  pathway  [6,12-15]. 
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Author  Summary 

Telomeres,  the  protective  elements  at  the  ends  of 
chromosomes,  need  to  be  maintained  for  cells  to 
proliferate  indefinitely.  In  many  human  cancers,  the 
telomeric  DNA  is  replenished  by  telomerase.  However,  a 
second  pathway  for  telomere  maintenance,  referred  to  as 
the  ALT  pathway,  has  increasingly  been  recognized  in 
human  cancers.  The  genetic  basis  for  activation  of  ALT  is 
not  known,  but  recent  data  have  implicated  a  chromatin 
remodeling  complex  (ATRX/DAXX)  and  the  histone  variant 
H3.3  as  players  in  the  repression  of  ALT.  We  have 
examined  a  large  panel  of  ALT  cell  lines  for  their  genetic 
and  cell  biological  features  and  found  that  loss  of  ATRX  is  a 
common  event  in  the  genesis  of  ALT  lines.  In  addition,  we 
document  that  ALT  cell  lines  frequently  undergo  chromo¬ 
somal  changes  and  are  impaired  in  their  ability  to  detect 
and  repair  damage  in  their  DNA.  These  hallmarks  of  ALT 
are  expected  to  facilitate  the  detection  of  ALT-type  tumors 
in  the  clinic  and  may  lead  to  ALT-specific  treatments. 


ALT  was  discovered  in  a  subset  of  immortalized  cell  lines  that 
emerge  at  low  frequency  from  human  cell  cultures  in  telomere 
crisis,  but  this  pathway  has  increasingly  been  identified  in  human 
cancer.  Approximately  10-15%  of  human  cancers  are  suspected  to 
use  the  ALT  pathway  based  on  their  lack  of  telomerase  activity  in 
combination  with  certain  hallmarks  of  ALT,  such  as  the  presence 
of  extrachromosomal  telomeric  circles,  APBs,  and/or  long  and 
heterogeneous  telomeres.  Based  on  these  criteria,  ALT  is  relatively 
common  in  many  sarcomas  (osteosarcoma  and  some  types  of  soft 
tissue  sarcomas),  certain  endocrine  tumors  (pancreatic  neuroen¬ 
docrine  tumors,  paraganglioma),  a  subset  of  nervous  system 
tumors  (e.g.  glioblastoma,  medulloblastoma,  oligodendroglioma, 
astrocytoma,  ganglioneuroblastoma),  bladder  small  cell  carcino¬ 
ma,  and  nonseminomatous  germ  cell  tumors  [16-18]. 

The  molecular  basis  for  the  activation  of  the  ALT  pathway  is 
not  known.  In  vitro,  the  frequency  of  conversion  to  ALT  is  low, 
typically  requiring  many  months  of  culturing  of  virally-trans- 
formed  (p53/Rb  deficient)  human  cells  that  have  entered  telomere 
crisis  [1 1].  The  low  frequency  of  conversion  to  ALT  suggests  that 
mutations  and/ or  infrequent  epigenetic  alterations  are  required  to 
unleash  this  pathway.  Indeed,  in  telomerase-positive  mammalian 
cells  with  fully  functional  telomeres,  telomeric  recombination  is 
stringently  repressed  by  the  TRF2,  Rapl,  and  POTl  components 
of  shelterin  as  well  as  by  the  Ku70/80  heterodimer  [6,19-23]. 
Mutations  in  shelterin  and  Ku  might  therefore  be  anticipated  in 
ALT,  a  notion  that  is  tested  here. 

Recently,  a  subset  of  pancreatic  neuroendocrine  tumors 
(PanNETs)  were  found  to  have  ultrabright  telomeric  FISH  signals, 
suggesting  the  presence  of  ALT-like  long  telomeres  [24,25].  These 
tumors  also  exhibited  alterations  in  ATRX  or  its  binding  partner 
DAXX.  Twenty-five  PanNETs  with  ultrabright  telomeric  foci  had 
ATRX/DAXX  mutations  or  lacked  nuclear  ATRX/DAXX  pro¬ 
tein.  In  contrast,  16  PanNETs  with  apparently  normal  ATRX/ 
DAXX  lacked  the  aberrant  telomeric  staining  patterns.  The  same 
correlation  between  inactivation  of  the  ATRX  pathway  and  the 
ALT-like  phenotype  was  observed  in  pediatric  glioblastoma  and 
several  other  cancers  [16,26]. 

ATRX  and  DAXX  act  together  in  a  replication-independent 
chromatin  assembly  pathway  that  deposits  the  histone  variant 
H3.3  at  telomeres  and  perhaps  at  other  G-rich  repetitive  elements 
[27-31].  The  function  of  ATRX/DAXX  and  histone  H3.3  at 
telomeres  is  not  yet  clear,  but  diminished  ATRX  function  has 


been  documented  to  increase  TTAGGG  repeat  containing 
telomeric  transcripts  (referred  to  as  TERRA  [32]),  reduce 
telomeric  loading  of  HP  1  a,  and  cause  modest  levels  of  telomere 
dysfunction  in  mouse  ES  (but  not  NIH3T3)  cells  as  gleaned  from 
the  localization  of  Y-H2AX  at  chromosome  ends  [29,31].  ATRX 
deficiency  also  leads  to  a  defect  in  sister  chromatid  cohesion  and 
aberrant  mitoses,  culminating  in  the  formation  of  micronuclei, 
lobulated  nuclei,  and  chromatin  bridges  [33-35]. 

Because  functional  assessment  of  ALT  is  not  possible  in  tumor 
specimens,  further  tests  for  the  association  of  ALT  with  ATRX/ 
DAXX  deficiency  requires  analysis  of  in  vitro  immortalized  cell 
lines.  Here  we  describe  the  results  of  a  comprehensive  effort  to 
characterize  the  genetic  alterations  and  associated  phenotypes  of 
22  human  ALT  cell  lines.  We  found  loss  of  ATRX  in  90%  of  in 
vitro  immortalized  ALT  lines,  suggesting  that  inactivation  of 
ATRX  is  a  major  step  in  generating  the  ALT  phenotype.  We  also 
document  that  ALT  cell  lines  are  characterized  by  extensive 
genomic  instability,  formation  of  micronuclei,  an  aberrant  G2/M 
checkpoint,  and  abnormal  DSB  repair  kinetics. 

Results 

Most  human  ALT  cell  lines  lack  ATRX 

We  assembled  a  panel  of  22  human  ALT  cell  lines  (Table  1), 
nineteen  of  which  originated  from  in  vitro  immortalization 
experiments  and  four  of  which  were  derived  from  human  tumors 
(three  osteosarcomas  and  one  lung  adenocarcinoma).  Many  of  the 
in  vitro  immortalized  lines  expressed  SV40-  or  HPV-derived 
oncoproteins.  Seven  cell  lines  originated  from  a  single  individual 
with  cystic  fibrosis  (JFCF-6  series);  an  additional  five  cell  lines  were 
established  from  a  Li-Fraumeni  patient  (IIICF  series).  We 
confirmed  the  absence  of  telomerase  activity  by  TRAP  assay 
(data  not  shown)  and  the  presence  of  extra-chromosomal  telomeric 
C-circles,  supporting  the  interpretation  that  these  were  ALT  cells 
[8,17]  (Table  1;  Figure  SI). 

The  ALT  cell  lines  and  non- ALT  controls  (hTERT-immortal- 
ized  SV40-transformed  BJ  fibroblasts  and  the  telomerase-positive 
HPV-E6/E7  expressing  HeLa  cervical  tumor  cell  line)  were  tested 
for  the  expression  of  ATRX  and  DjAXX  (Figure  lA-lC  and  data 
not  shown).  Repeated  immunoblots  on  independently  harvested 
whole  cell  extracts  indicated  that  among  the  22  ALT  cell  lines, 
only  six  (G292,  IIICF-E6E7/C4,  IVG-BF.LXN,  LFS-05F-24,  SK- 
LU-1,  and  MeT-4A),  contained  detectable  full  length  ATRX 
protein  and  one  of  these  (IVG-BF.LXN)  appeared  to  lack  DAXX. 
Of  the  six  cell  lines  expressing  ATRX,  three  (G292,  IIICF-E6E7/ 
C4,  and  IVG-BF.LXN)  showed  very  low  immunofluorescence 
signals  for  ATRX  and/or  DAXX  (Figure  IB),  whereas  the  other 
three  lines  displayed  the  punctate  nuclear  staining  pattern  for 
ATRX  and  DTkXX  expected  for  this  PML  body  component 
[36,37].  Thus,  19  out  of  22  ALT  lines  had  an  alteration  in  the 
expression  of  ATRX  and/or  DAXX.  As  each  of  the  seven  JFCF- 
6-derived  ALT  lines  lacked  detectable  ATRX,  we  considered  the 
possibility  that  the  parental  JFCF-6  culture  might  have  been 
deficient  in  ATRX  expression.  However,  a  non-ALT  JFCF-6  line 
contained  ATRX  (Figure  1 C),  indicating  that  the  lack  of  ATRX  in 
the  other  JFCF-6  lines  was  potentially  correlated  with  the  ALT 
phenotype. 

All  ALT  cell  lines  were  tested  for  genetic  changes  in  ATRX, 
DAXX,  and  H3.3  {H3F3A)  using  primers  designed  to  PCR-amplify 
all  exons  from  these  genes  for  Sanger  sequencing  (Table  1;  Table 
SI).  Failure  to  amplify  multiple  consecutive  exons  was  interpreted 
as  a  deletion  event  only  when  flanking  exons  did  amplify.  By  that 
criterion,  6  cell  lines  gFCF-6/T.  ID,  JFCF-6/T.  IL,  JFCF-6/T.  IR, 
JFCF-6/T.1J/5H,  JFCF-6/T.1J/1-3C,  and  IIICF-T/B3)  har- 
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bored  large  deletions  in  the  ATRX  gene  that  ranged  from  4  to  26 
exons  in  size  (Table  1;  Table  Si). 

Sequencing  of  genes  that  did  not  harbor  detectable  deletions 
revealed  no  potentially  inactivating  genetic  changes.  Nine  of  those 
cell  lines  lacked  detectable  ATRX  protein  suggesting  that  genetic 
alterations  not  detectable  by  Sanger  sequencing,  such  as  translo¬ 
cations,  promoter  or  splicing  mutations,  or  epigenetic  changes  may 
underlie  their  lack  of  expression.  No  genetic  alterations  were 
detected  in  DAXX  or  the  H3F3A  gene  encoding  H3.3  (Table  1). 

These  data  provide  strong  evidence  that  loss  of  ATRX 
expression  is  involved  in  either  the  initiation  or  maintenance  of 
the  ALT  pathway  in  human  cells.  To  test  whether  deficiency  in 
ATRX  was  sufficient  to  induce  the  telomere-telomere  recombi¬ 
nation  typical  of  ALT  cells,  we  used  an  shRNA  to  deplete  ATRX 
from  HeLa  cells  and  used  Chromosome-Orientation  FISH  (CO¬ 
FISH)  to  measure  sequence  exchanges  between  sister  telomeres 
(Telomere  Sister  Chromatid  Exchanges  or  T-SCEs).  Despite 
repression  of  ATRX  by  more  than  95%,  the  frequency  of  T-SCEs 
was  not  significantly  altered  (Figure  ID  and  IE).  In  addition, 
repression  of  ATRX  or  DAXX  in  SV40-transformed  BJ 
fibroblasts  failed  to  induce  escape  from  crisis  (Figure  IF  and  IG; 
Figure  S2),  whereas  expression  of  hTERT  readily  immortalized 
the  cells  (Figure  IG;  Figure  S2).  After  a  prolonged  period  in  crisis 
(6  weeks),  immortal  clones  emerged  from  the  culture  in  which 
ATRX  expression  was  suppressed  with  ATRX  sh590.  However,  a 
parallel  culture  in  which  ATRX  was  not  suppressed  efficiently 
(ATRX  sh589)  also  yielded  an  immortal  subpopulation  (Figure  IF 
and  1 G)  and  in  both  cases,  the  immortal  cells  expressed  telomerase 
and  lacked  telomeric  C-eircles  (Figure  S2A  and  S2B),  indicating 
that  they  had  not  activated  AET.  These  observations  argue  that 
this  level  of  ATRX/DAXX  suppression  is  not  suffieient  to  activate 
the  AET  pathway. 

No  overt  changes  in  shelterin  or  299  potentially  ALT 
relevant  genes 

The  lack  of  induction  of  AET  by  shRNA-mediated  suppression 
of  ATRX  and  DAXX  prompted  us  to  inspect  the  AET  cell  line 
panel  for  changes  in  the  status  of  genes,  pathways,  and  proteins 
with  possible  relevance  to  the  activation  of  AET.  Because  HDR  is 
normally  repressed  at  telomeres  by  shelterin,  we  examined  the 
expression  of  all  six  shelterin  components  in  the  AET  lines  by 
repeated  immunoblotting  (Figure  S3).  Using  this  approach  on  all 
22  ALT  lines  failed  to  reveal  a  consistent  change  in  the  abundance 
or  appearance  of  the  six  shelterin  proteins.  In  addition,  exome 
sequencing  (see  below)  failed  to  reveal  mutations  in  the  genes 
encoding  TRFl,  TRF2,  Rapl,  POTl,  TIN2,  and  TPPl. 

A  second  entity  associated  with  telomeres  is  the  telomeric  RNA, 
TERRA  [32].  As  mouse  ATRX  deficient  cells  were  reported  to 
have  elevated  TERRA  levels  and  selected  AET  cells  were 
previously  found  to  have  higher  levels  of  TERRA  [31,38],  we 
examined  the  expression  of  this  RNA  in  the  AET  lines. 
Quantitative  Northern  analysis  confirmed  that  the  AET  lines 
generally  have  higher  TERRA  levels  than  telomerase-positive 
cells,  but  some  AET  lines  (e.g.  JFCF-6/T.  ID  and  IIICF-E6E7/ 
C4)  approximated  the  levels  of  TERRA  found  in  BJ/hTERT/ 
SV40  and  HeEa  cells  (Figure  S4).  There  was  no  clear  correlation 
between  increased  TERRA  and  ATRX  status  since  the  cell  lines 
with  apparendy  normal  ATRX/DAXX  (EFS-05F-24,  SK-EU-1, 
and  MeT-4A)  had  TERRA  levels  comparable  to,  or  higher  than, 
cell  lines  lacking  ATRX  (e.g.  JFCF-6/T.1D,  JFCF-6/T.1R,  and 
JFCF-6/T.1E)  (Figure  S4). 

In  addition  to  shelterin  and  TERRA,  we  analyzed  299  genes 
with  potential  relevance  to  telomeres  and  the  AET  pathway.  This 
gene  set  included  genes  encoding  proteins  involved  in  telomere 
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Figure  1.  Deficiency  in  ATRX/DAXX  correlates  with  ALT.  A,  Immunoblot  for  ATRX  and  DAXX  in  the  indicated  ALT  lines.  ATRX  was  detected 
with  A30T045  (Bethyl  Labs).  DAXX  was  detected  with  A301-353A  (Bethyl  Labs).  Asterisk:  cell  lines  with  deletions  in  the  ATRX  gene  (see  Table  1  and 
Table  S2).  HeLa  cells  (blue)  were  used  as  a  positive  control.  B,  IF  for  ATRX  and  DAXX  in  the  indicated  cell  lines.  HeLa  cells  (blue)  are  used  as  a  positive 
control.  C,  Immunoblot  for  ATRX  in  additional  non-ALT,  positive  control  cells  (blue)  relative  to  the  indicated  ALT  lines.  D,  Immunoblot  showing  >90% 
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reduction  of  ATRX  protein  in  HeLa  cells  expressing  ATRX  shRNA592.  E,  T-SCEs,  a  measure  for  telomeric  recombination,  assessed  by  Chromosome 
Orientation  (CO)-EISEI  on  metaphases  harvested  from  EleLa  cells  expressing  vector  or  ATRX  shRNA592.  The  average  percentage  of  T-SCEs  from  four 
independent  experiments  is  shown,  with  a  p  value  derived  from  a  two-tailed,  unpaired  f  test.  E,  Immunoblot  showing  >90%  reduction  of  ATRX 
protein  in  cells  expressing  ATRX  shRNA590  but  not  with  shRNA589.  G,  Growth  curves  showing  immortalization  of  SV40-transformed  BJ  fibroblasts 
infected  with  an  hTERT  expressing  retrovirus  but  no  immortalization  after  infection  with  the  effective  (sh590)  and  ineffective  (sh589)  shRNAs  to  ATRX. 
doi:10.1371/journal.pgen.1002772.g001 


maintenance,  DNA  damage  repair  and  signaling,  chromosome 
duplication  and  segregation,  and  cell  cycle  regulation  (Table  S2). 
The  coding  regions  of  these  genes  were  analyzed  by  massively 
parallel  454  sequencing  on  genomic  DNA  derived  from  1 4  of  the 
22  ALT  lines.  We  focused  on  ALT  lines  generated  in  vitro  as 
opposed  to  tumor-derived  ALT  lines  to  exclude  genetic  changes 
selected  for  during  tumorigenesis.  The  JFCF-6  and  IIICF  ALT 
lines  are  derived  from  two  individuals,  allowing  identification  of 
SNPs  specific  to  the  individual  as  well  as  potential  ALT-correlated 
changes.  After  excluding  known  SNPs,  the  remaining  sequence 
data  were  filtered  to  exclude  apparent  SNPs  present  in  less  than 
25%  of  the  reads.  The  rare  SNPs  that  eluded  the  threshold  filter 
affected  20  genes  (Table  S3).  Within  the  JFCF-6  and  IIICF  cell 
line  sets,  a  nucleotide  change  in  multiple  members  of  either  set  was 
interpreted  as  being  present  prior  to  the  acquisition  of  the  ALT 
phenotype,  and  thus  not  relevant.  The  variable  occurrence  of 
these  SNPs  in  the  multiple  ALT  lines  derived  from  one  individual 
(e.g.  only  two  of  the  five  IIICF  derived  lines  have  the  BLM 
K323R)  is  likely  due  to  chromosome  losses  occurring  during  or 
after  telomere  crisis.  The  florid  chromosome  instability  exhibited 
by  these  cell  lines  (see  below)  makes  this  an  appealing  and 
parsimonious  explanation.  Some  of  these  SNPs  were  also  found  in 
other  cell  lines  (e.g.  MDCl  C1559G  in  LSF-05F-24  and  IVG- 
BF.LXSN)  and  deemed  irrelevant  based  on  their  presence  in  the 
IIICF  and/or  JFCF-6  lines  prior  to  the  induction  of  ALT.  Based 
on  these  criteria,  no  ALT  specific  change  was  observed  in  the  five 
IIICF-derived  ALT  lines  or  in  IVG-BF.LXN.  Three  of  the  JFCF-6 
derived  ALT  lines  had  a  specific  change  that  could  not  be 
attributed  to  a  pre-existing  SNP  (A436S  in  Suv420H2,  S553F/ 
S470F  in  PRMT5,  and  D2579V  in  Ki67;  Table  S3).  Since  these 
genes  were  not  affected  in  the  other  13  ALT  lines  analyzed,  they 
are  unlikely  to  be  critical  to  ALT.  Similarly,  LSF-05F-24 
contained  changes  in  five  genes  that  are  not  affected  in  the  other 
13  cell  lines.  Since  none  of  these  alterations  were  found  in  a 
substantive  fraction  of  the  ALT  cell  lines,  it  is  unlikely  that  these 
genes  represent  clear  contributors  to  the  ALT  phenotype. 

Extensive  genome  rearrangements  in  ALT 

As  ALT-relevant  genetic  alterations  might  be  identified  from 
recurrent  regions  of  gene  copy  number  changes,  we  analyzed  the 
ALT  cell  line  panel  using  Affymetrix  6.0  single  nucleotide 
polymorphism  (SNP)  arrays  to  identify  regions  of  recurrent  copy 
number  alteration.  When  we  examined  either  individual  cells 
(Figure  2A  and  2B)  or  the  entire  set  as  a  group  (GISTIC  analysis. 
Figure  S5),  we  found  that  the  number  and  complexity  of  these 
gains  and  losses  were  at  least  as  extensive  as  the  alterations  in 
human  cancers  with  genome  instability  [39].  When  we  considered 
all  of  the  ALT  cell  lines  as  a  group,  we  identified  18  regions  of 
significant  focal  copy  number  gain,  and  many  more  regions  of 
copy  number  loss  (Figure  S5). 

The  genome  rearrangements  in  the  ALT  lines  are  a  likely 
consequence  of  the  telomere  dysfunction  these  cells  experienced 
prior  to  their  immortalization.  Consistent  with  this  idea,  SKY 
analysis  of  five  ALT  lines  JFCF-6/T.1R,  JFCF-6/T.1Q,  IVG- 
BF.LXN,  IIICF/c,  and  LFS-05F-24)  showed  frequent  non¬ 
reciprocal  translocations,  deletions,  complex  rearrangements, 
and  hyper-triploid  chromosome  numbers  (Figure  3),  all  of  which 


can  be  caused  by  telomere  dysfunction  [40,41].  Previous 
cytogenetic  analysis  of  two  ALT  cell  lines  (Saos-2  and  ZK-58) 
also  indicated  frequent  translocations,  deletions,  and  complex 
rearrangements  [42]. 

In  cells  that  have  been  immortalized  by  telomerase  activation, 
the  telomere-driven  genome  instability  is  largely  dampened, 
resulting  in  a  rearranged  but  now  relatively  stable  karyotype 
[43,44].  We  asked  whether  ALT  similarly  stabilizes  the  genome  by 
examining  two  generations  of  clonal  descendants  of  the  JFCF-6/ 
T.  IR  ALT  line.  Although  the  clones  exhibited  features  in  common 
with  the  parental  cell  line,  each  of  the  clones  showed  new  regions 
of  copy  number  gain  or  loss  (Figure  2B).  For  instance,  a  new  gain 
of  a  segment  on  chromosome  2q  was  observed  in  one  single  cell 
clone  of  the  parental  line  (and  all  the  clones  derived  from  it),  and 
in  another  clone  (and  its  derivatives)  an  overlapping  segment 
showed  a  new  deletion.  Many  other  changes,  including  examples 
on  chromosomes  6,  12,  14,  15,  16  and  21,  showed  obvious  new 
gains  or  losses  as  compared  to  the  parental  clone,  with 
independent  clones  showing  different  alterations.  In  contrast,  we 
found  only  rare  clonal  alterations  in  single  cell  clones  derived  from 
the  telomerase-expressing  PICT  116  cell  line.  These  observations 
suggest  that  the  JFCF-6/T.  IR  and  perhaps  other  ALT  lines  have 
ongoing  genome  instability  as  well  as  highly  rearranged  genomes. 

We  considered  that  ongoing  genome  instability  in  the  ALT  lines 
could  be  due  to  defects  in  mitosis  and/or  dicentric  chromosomes 
resulting  from  telomere  dysfunction  and  ATRX  deficiency. 
Defects  in  mitosis  and  dicentric  chromosomes  can  give  rise  to 
lagging  chromosomes  that  can  form  micronuclei,  which  were 
recently  shown  to  engender  massive  genomic  alterations  [45].  We 
therefore  examined  the  ALT  lines  for  the  spontaneous  occurrence 
of  micronuclei.  Fifteen  of  the  22  ALT  lines  showed  a  high 
frequency  (10-30%)  of  cells  with  micronuclei  (Table  1;  Figure  4A 
and  4B).  As  expected,  SV40-transformed  hTERT-expressing  BJ 
fibroblasts  and  HeLa  cells  showed  the  low  level  of  spontaneous 
micronucleation  (<8%)  previously  noted  in  other  transformed 
human  cells  [46].  Micronucleation  frequencies  of  up  to  15%  were 
previously  only  observed  in  cells  experiencing  high  levels  of  DNA 
damage  ensuing  from  yirradiation,  in  genetic  contexts  associated 
with  high  levels  of  spontaneous  DNA  damage,  and  in  genomically 
unstable  tumor  cell  lines  (e.g.  [47,48];  Neill  Ganem  and  David 
Pellman,  pers.  comm.).  Thus,  many  ALT  cell  lines  show  an 
unusually  high  frequency  of  micronucleation,  indicative  of 
ongoing  genome  instability.  The  micronucleation  phenotype  of 
ALT  lines  is  likely  to  be  in  part  due  to  the  absence  of  ATRX  since 
depletion  of  ATRX  with  two  shRNAs  induced  the  formation  of 
micronuclei  in  BJ/hTERT/SV40  and  HeEa  cells  (Figure  4C  and 
4D),  consistent  with  a  previous  report  on  the  induction  of 
lobulated  nuclei  after  ATRX  depletion  [35]. 

G2/N\  checkpoint  function  and  DSB  repair  kinetics  in  ALT 
lines 

The  signs  of  ongoing  genome  instability  in  ALT  cells  led  us  to 
examine  the  functionality  of  DNA  damage  checkpoints  and  DSB 
repair.  As  the  ALT  lines  were  expected  to  have  an  impaired  G 1  /  S 
DNA  damage  checkpoint  due  to  their  p53  deficiency,  we  focused 
our  efforts  on  the  G2/M  checkpoint.  Entry  into  mitosis  is  blocked 
upon  activation  of  the  ATM  and/or  ATR  signaling  pathways. 


PLoS  Genetics  |  www.plosgenetics.org 


6 


July  2012  I  Volume  8  |  Issue  7  |  el  002772 


ATRX  Loss  and  Genome  Instability  in  ALT 


Cell  line 
IIICF-T/A6 
IIICF-T/B3 
IIICF/a2 
IIICF/c 
IIICF-E6E7/C4 
JFCF-6/T.1J/1-3C 
JFCF-6/T.1J/5H 
JFCF-6/T.1Q 
JFCF-6/T.1M 
JFCF-6/T.1R 
JFCF-6/T.1L 
JFCF-6/T.1D 
G292 
LSF-05F-24 
SK-LU-1 
SUSM-1 
U-2  OS 
GM847 
MeT-4A 
WI38-VA12/2RA 
IVG-BF-LXSN 
SaOS-2 


1 


g  Chromosome 

Cell  line  ^1,2,3,  4  ,  5  ,  6  ,  7  ,  8  ,  9  ,  10 ,  11  , 12  ,13  ,14 ,15,16,^^,18'|^4  ^  X 

f-  HCT116  p 
HCT116/C1 
HCT116/C4 
HCT116/C5 
-JFCF-6/T.1R  p 
JFCF-6/T.1R/C3 
JFCF-6/T.1R/C7 
I-  LjFCF-6/T.1R/c8 


lUI  i 


Cell  line 
'-p  JFCF-6/T.1R/C8 
-JFCF-6/T.1R/c8c2 
-JFCF-6/T.1R/c8c5 
LjFCF-6/T.1R/c8c1 
'-p  JFCF-6/T.1R/C7 
■JFCF-6/T,1R/c7c1 
-JFCF-6/T,1R/c7c2 
'-JFCF-6/T.1R/c7c3 


Chromosome 


17  19  21 

18  20  22 

I  1  I  I  I  I  I 


X 

T 

U 


Y 

j  L 


Figure  2.  Copy  number  analysis  showing  extensive  genome  rearrangements  in  ALT  iines.  SNP  array  copy  number  results  are  shown  for 
(A)  the  ALT  cell  lines  and  (B)  the  1st  and  2nd  generation  subclones  derived  from  JFCF-6/T.1R  (indicated  by  the  brackets)  and  subclones  derived  from 
the  telomerase-positive  FICT1 16  control.  Segmented  copy  number  data  is  shown  for  each  chromosome,  by  genomic  position  in  columns  and  by  cell 
line  in  rows.  The  color  scale  ranges  from  red  (amplification;  log2  copy  number  ratio  of  1.5)  through  white  (neutral;  2  copies  in  diploid  lines,  log2  ratio 
of  0)  to  blue  (deletion;  log2  ratio  of  -1.5). 
doi:1 0.1 371/journal.pgen.1 002772.g002 


ATM  signaling  is  primarily  responsible  for  the  initiation  of  the 
checkpoint,  whereas  ATR  signaling  ensures  the  maintenance  of 
the  arrest  [49].  To  determine  the  efficiency  of  initiation  and 
maintenance  of  the  G2/M  checkpoint,  subconfluent  cells  were 


subjected  to  ionizing  radiation,  and  the  mitotic  index  was 
determined  by  flow  cytometric  quantitation  of  histone  H3  (SIO) 
phosphorylation.  One  hour  after  treatment  with  10  Gy  ionizing 
radiation  (IR),  the  mitotic  index  of  control  hTERT-immortalized 
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Figure  3.  Abnormal  karyotypes  in  ALT  lines.  Representative  metaphases  from  5  cell  lines  showing  subtetraploid  karyotype  with  high  number  of 
rearranged  chromosomes  (up  to  60%).  Structural  rearrangements  are  labeled  as  following:  del(Z)-  deletion  of  chromosome  Z;  der(Z)  -  multiple 
aberrations  within  single  chromosome  Z;  chromosome,  denoted  with  two  or  more  numbers  indicate  rearrangement  involving  two  or  more  partners. 
doi:10.1371/journal.pgen.1002772.g003 
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Figure  4.  Frequent  micronucleation  in  ALT  celis  may  be  attributabie  to  ioss  of  ATRX.  A,  Examples  of  micronuclei  (arrowheads)  in  two  ALT 

cell  lines.  Blue:  DAPI  stain  for  DNA.  B,  Graph  showing  micronucleation  frequencies  in  the  ALT  cell  lines  and  two  telomerase  positive  controls  (blue). 
Values  are  means  ±SD  from  three  experiments  (>100  nuclei  each).  Red:  cell  lines  with  >10%  micronucleation  frequency  (dashed  line).  HeLa  cells 
were  analyzed  once.  C,  Examples  of  micronuclei  (arrowheads)  in  HeLa  cells  expressing  two  independent  ATRX  shRNAs  (see  figure  1  for  immunoblots). 
D,  Graph  showing  micronucleation  frequencies  in  two  telomerase  positive  cell  lines  infected  with  vector  or  the  indicated  ATRX  shRNAs.  Values  are 
means  ±SEIVI  from  two  experiments  (200  nuclei  each). 
doi:10.1371/journal.pgen.1002772.g004 


Retinal  Pigment  Epithelial  cells  (RPE/hTERT)  was  reduced 
approximately  ten-fold  (from  2.7%  to  0.3%),  whereas  reduction  in 
the  mitotic  index  of  ATM-deficient  GM5849  cells  was  less  than 


two-fold  (from  4.1%  to  2.3%)  (Table  1;  Figure  5A).  Of  fourteen 
ALT  lines  that  had  a  sufficiently  high  mitotic  index  to  yield 
interpretable  results,  half  showed  a  severe  defect  in  G2/M 
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Figure  5.  Many  ALT  lines  have  a  defect  in  the  G2/M  checkpoint  initiation  and/or  maintenance.  A,  Bar  graph  depicting  the  results  of  an 
assay  for  the  initiation  of  the  G2/M  checkpoint  at  1  hr  after  irradiation  of  the  indicated  cell  lines.  RPE/hTERT  is  a  positive  control  (in  blue).  Mean  ±  SD 
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for  triplicate  experiments  are  shown.  B,  Bar  graph  depicting  the  results  of  an  assay  for  the  maintenance  of  the  G2/IVI  checkpoint  at  1 6  hr  after  IR.  BJ/ 
hTERT/SV40  is  a  positive  control  (in  blue).  Mean  ±  SD  for  triplicate  experiments  are  shown.  See  Table  1  for  summary.  Cell  lines  in  red  have  a  defect  in 
the  G2/M  checkpoint. 
doi:1 0.1 371/journal.pgen.1 002772.g005 


checkpoint  initiation  (Figure  5A;  Table  1).  In  six  ALT  lines,  the 
checkpoint  defect  was  more  severe  than  that  of  ATM-deficient 
cells.  Qualitatively  similar  results  were  observed  in  response  to  low 
dose  (1  Gy)  IR  treatment,  with  the  same  cell  lines  exhibiting 
defects  as  at  the  high  dose  (data  not  shown).  Similarly,  10  of  11 
ALT  lines  tested  displayed  a  deficiency  in  the  maintenance  of  the 


G2/M  checkpoint,  monitored  at  16  hr  after  y-irradiation  with  4.5 
Gy  (Table  1;  Figure  5B).  Some  of  the  ALT  lines  were  proficient  in 
the  initiation  of  the  checkpoint  but  failed  to  maintain  it  (e.g., 
IIICF/c).  Collectively,  we  observed  impaired  G2/M  checkpoint 
function  in  11  out  of  1 7  ALT  lines  tested.  This  defect  in  the  G2  /M 
checkpoint  is  not  seen  in  established  non-ALT  human  cell  lines. 
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Figure  6.  Deficiency  in  DSB  repair  in  ALT  iines.  A,  Example  of  assay  for  53BP1  foci  in  the  positive  control  (BJ/hTERT/SV40,  blue)  (top)  and  an  ALT 
line  (bottom).  Cells  were  treated  as  indicated  above  the  panels  and  stained  for  DNA  with  DAPI  (blue)  in  conjunction  with  IE  for  53BP1  (green).  Note 
the  higher  level  of  53BP1  foci  in  non-irradiated  ALT  cells  and  incomplete  DSB  repair  at  24  hr  after  IR.  B,  Example  of  PEGE-assay  for  DSB  repair  in  the 
indicated  cell  lines.  Cells  were  treated  with  100  Gy  IR  24  hr  before  harvesting  and  residual  DSBs  were  evaluated  based  on  the  DNA  fragments  released 
from  agarose  embedded  cells  into  a  PEGE  gel.  Plugs  were  treated  with  600  Gy  to  fragment  the  DNA  and  run  in  parallel  to  serve  as  a  control  for 
loading.  DNA  was  stained  with  ethidium  bromide  after  electrophoresis.  HT1080  (fibrosarcoma  cell  line)  was  used  as  a  positive  control  (in  blue). 
doi:1 0.1 371/Joumal.pgen.1 002772.g006 
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For  example,  normal  G2/M  arrest  following  y-irradiation  is 
observed  in  HT1080,  HCT116,  IMR90  gHJP,  unpubl.  data), 
293T  [50],  MCF7,  and  genetically  complemented  HCC1937  cells 
[51].  The  decrement  in  the  G2/M  checkpoint  did  not  reflect  a 
global  impairment  of  DNA  damage  response  (DDR)  signaling  in 
the  ALT  cells,  since  phosphorylation  of  CFIK2,  an  outcome  of 
ATM  signaling,  was  readily  apparent  upon  IR  treatment  in  all  cell 
lines,  but  undetectable  in  GM5849  (ATM-deficient)  cells  (Figure 
S6).  Indeed,  CHK2  phosphorylation  was  evident  in  several  of  the 
ALT  cell  lines  prior  to  irradiation,  consistent  with  the  finding  that 
the  ALT  phenotype  is  associated  with  chronic  genotoxic  stress 
(Figure  S6). 

In  addition  to  an  abnormal  G2/M  checkpoint,  many  of  the 
ALT  cell  lines  showed  unusually  slow  DSB  repair  kinetics.  To 
assay  this  aspect  of  the  DDR,  cells  were  y-irradiated  with  0.5  Gy 
and  the  fraction  of  cells  with  10  or  more  53BP1  foci  was  scored 
after  1  and  24  hr  as  well  as  in  non-irradiated  controls  (Figure  6A; 
Figure  S7;  Table  1).  As  expected,  non-ALT  cell  lines  (HeLa  and 
hTERT  immortalized  SV40-transformed  BJ)  showed  a  nearly 
complete  disappearance  of  the  IR-induced  DNA  damage  foci  at 
24  hr.  In  contrast,  1 1  out  of  12  ALT  cell  lines  retained  a  significant 
fraction  of  the  induced  53BP1  DNA  damage  foci  after  a  day 
(Table  1;  Figure  6A).  This  result  was  corroborated  with  a  pulsed 
field  gel  electrophoresis  (PFGE)-based  gel  assay  for  the  disappear¬ 
ance  of  IR-induced  DNA  fragmentation  after  24  hr  (Figure  6B; 
Table  1).  We  also  noted  that  approximately  half  of  the  ALT  cell 
lines  had  a  high  basal  level  of  cells  containing  10  or  more  53BP1 
foci,  further  confirming  that  ALT  cells  experience  ongoing 
genome  damage  (Table  1).  These  spontaneous  53BP1  foci  could 
be  the  result  of  both  genome-wide  DNA  breaks  and  the  presence 
of  dysfunctional  telomeres  [52].  Their  deficiency  in  the  G2/M 
checkpoint,  combined  with  the  absence  of  a  functional  p53 
pathway  likely  explains  why  these  cells  proliferate  despite  severe 
genome  damage. 

Discussion 

Here  we  report  that  ATRX  is  either  undetectable  or  severely 
depleted  from  PML  bodies  in  ~90%  of  human  ALT  cell  lines 
tested,  including  ALT  lines  derived  in  vitro.  This  result  establishes 
a  strong  correlation  between  the  initiation  or  maintenance  of  ALT 
and  deficiency  in  the  ATRX/DAXX  pathway.  In  a  subset  of 
cases,  the  loss  of  ATRX  is  due  to  large  deletions  within  the  ATRX 
gene,  further  underscoring  the  relevance  of  ATRX  status  to  ALT. 
We  expect  that  alterations  in  ATRX/DAXX  in  cells  and  tumors 
will  be  a  useful  indicator  of  the  presence  of  the  ALT  pathway. 
However,  manipulating  ATRX/DAXX  expression  failed  to 
unleash  ALT,  suggesting  that  deficiency  of  these  proteins  is  not 
sufficient  for  the  ALT  phenotype  and  pointing  to  the  need  to 
identify  the  cooperating  (epi-)  genetic  changes.  In  addition,  a  key 
issue  for  further  investigation  is  the  mechanism  by  which  the  loss 
of  ATRX/DAXX  and  concomitant  lack  of  deposition  of  H3.3 
into  telomeric  chromatin  might  allow  or  promote  telomere 
recombination.  We  note  that  it  remains  possible  that  the 
ATRX/DAXX  pathway  has  an  indirect  effect,  for  instance  by 
repressing  genes  that  promote  telomere  recombination. 

Our  data  have  uncovered  a  hitherto  unappreciated  level  of 
genome  instability  in  ALT  cells.  All  ALT  lines,  including  cell  lines 
that  arose  from  in  vitro  immortalization  experiments,  harbor 
scrambled  genomes  and  have  signs  of  ongoing  genome  instability. 
The  ALT  lines  have  frequent  micronuclei,  a  high  basal  level  of 
DNA  damage  foci,  and  elevated  checkpoint  signaling  in  absence  of 
exogenous  damage.  Furthermore,  in  one  ALT  cell  line,  repeated 
subcloning  established  frequent  copy  number  alterations.  These 


features  of  ALT  raise  the  possibility  that  tumors  employing  this 
pathway  will  have  unique  vulnerabilities  that  could  offer  clinical 
utility,  such  as  has  been  shown  for  cells  lacking  BRCAl  that  are 
sensitive  to  PART  inhibition  [53].  This  possibility  underscores  the 
importance  of  establishing  a  signature  for  the  ALT  phenotype. 

In  this  context,  it  is  notable  that  G2/M  checkpoint  deficiency 
emerged  as  a  prominent  attribute  of  many  ALT  cells  in  this  study. 
It  is  likely  that  the  diminished  G2/M  checkpoint  response  allows 
these  cells  to  proliferate,  despite  a  considerable  burden  of  DNA 
damage  both  at  telomeres  and  elsewhere  in  the  genome.  G2/M 
checkpoint  deficiency  is  an  example  of  a  potential  vulnerability 
unique  to  ALT  cells.  Indeed,  G2/M  checkpoint  inhibitors  to 
enhance  efficacy  of  clastogenic  therapies  are  currently  being 
developed  and  evaluated  in  clinical  trials  [54]. 

The  altered  G2/M  response  could  not  be  ascribed  to  mutations 
in  known  components  of  the  ATM  or  ATR  signaling  pathways. 
Furthermore,  DNA  damage  signaling  was  demonstrably  intact 
since  CHK2  phosphorylation  could  be  induced  by  exogenous 
DNA  damage.  The  implication  of  ATRX  and  DAXX  in  the  ALT 
phenotype  may  have  some  relevance  to  the  checkpoint  defects 
observed,  as  the  involvement  of  chromatin  remodelers  in  the  DDR 
becomes  increasingly  clear  [55].  Nevertheless,  identification  of  the 
mechanisms  underlying  G2/M  checkpoint  dysfunction  in  ALT 
cells  is  warranted  given  the  data  presented  here. 

Our  results  would  appear  to  exclude  the  possibility  that  the 
ALT  cells  arise  from  a  single  dominant  mutation  in  addition  to 
p53/Rb  loss.  Loss  of  ATRX/DAXX  function  appears  to  be 
required  for  ALT  but  it  seems  likely  that  a  defect  in  the  G2/M 
checkpoint  is  also  needed  in  order  for  the  cells  to  proliferate,  given 
their  high  level  of  spontaneous  DNA  damage.  The  view  that 
multiple  steps  are  required  to  allow  ALT-mediated  immortaliza¬ 
tion  is  consistent  with  the  low  frequency  by  which  it  arises,  even  in 
SV40-  or  HPV-transformed  cells. 

Materials  and  Methods 

Cell  lines 

ATL  cell  lines  are  described  in  Table  1.  CellBank  Australia 
provided  cell  line  quality  control.  All  cell  lines  were  verified  by  1 6- 
locus  STR  profiling  and  confirmed  to  be  free  of  Mycoplasma 
species,  or  were  analyzed  within  1 0  population  doublings  of  being 
obtained  from  ATCC.  Non-ALT  cell  lines  were  obtained  from  the 
ATCC.  The  HeLa  cell  line  (subclone  HeLal.3)  was  described 
previously  [56]. 

Sequence  analysis  of  ATRX,  DAXX,  and  H3.3  genes 

Primers  were  synthesized  by  Invitrogen  (San  Diego,  CA)  to 
cover  all  exons  oi  ATRX,  DAXX,  and  H3F3A1  (H3.3).  PGR  was 
performed  in  5  pi  reactions  containing  1  x  PGR  Buffer  (67  mM 
Tris-HCl,  pH  8.8,  6.7  mM  MgC12,  16.6  mM  NH4S04,  10  mM 
2-mercaptoethanol),  1  mM  dNTPs  (Invitrogen,  San  Diego,  CA), 
1  pM  forward  and  1  pM  reverse  primers,  6%  DMSO,  2  mM 
ATP,0.25  U  Platinum  Taq  (Invitrogen,  San  Diego,  CA)  and  3  ng 
DNA.  Reactions  were  carried  out  in  384-well  ABI  9700 
thermocyclers  (Applied  Biosystems,  Foster  City,  CA)  using  a 
touchdown  PCR  protocol:  1  cycle  of  96°C  for  2  min;  3  cycles  of 
96°C  for  10  sec,  64°C  for  10  sec,  70°C  for  30  sec;  3  cycles  of  96°C 
for  10  sec,  61°C  for  10  sec,  70°C  for  30  sec;  3  cycles  of  96°C  for 
10  sec,  58°C  for  10  sec,  70°C  for  30  sec;  41  cycles  of  96°C  for 
10  sec,  57°C  for  10  sec,  70°C  for  30  sec;  Icycle  of  70°C  for  5  min. 
The  PCR  products  were  evaluated  for  presence  and  size  by 
electrophoresis  on  2%  agarose  gels  and  sequenced  as  described 
(refs.  25  and  57  [25,57]).  All  failed  PCRs  were  repeated  at  least 
two  more  times. 
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Targeted  sequencing  of  candidate  ALT  genes 

Coding  regions  of  299  genes  were  amplified  from  genomic 
DNA  of  15  samples.  Following  PCR  amplification,  products  were 
normalized  and  pooled  prior  to  emulsion  PCR  amplification  of 
single  stranded  DNA  and  454  pyrosequencing.  Sequencing  was 
performed  using  a  FLX  Titanium  machine  at  Agencourt 
Bioscience  Corp.  (Beverly  MA),  to  a  depth  of  20  x  coverage  per 
sample  per  amplicon  (attempted).  Assembly  and  mapping  of  reads 
was  performed  using  AVA  software. 

Immunoblotting 

Immunoblotting  was  performed  using  standard  procedures  on 
whole  cell  extracts.  ATRX  was  detected  primarily  with  A30 1-045 
(Bethyl  Labs)  or  with  FIPA001906  (Sigma  Aldrich).  DAXX  was 
detected  with  A301-353A  (Bethyl  Labs).  Immunoblotting  for 
CHK2-P  was  done  with  CHK2-T69  (Cell  Signaling)  and  total 
CHK2  was  detected  with  clone  7  Ab  (Millipore).  Shelterin 
components  were  analyzed  as  described  previously  [56]. 

Immunofluorescence 

Cells  grown  on  glass  coverslips  were  fixed  with  3%  parafor¬ 
maldehyde  for  10  min  at  room  temperature,  washed  three  times 
with  PBS,  and  permeabilized  with  0.5%  Triton  X-100  in  PBS  for 
10  min  on  ice.  Cells  were  then  washed  four  times  with  PBS  and 
blocked  for  30  min  at  room  temperature  with  PBG  (0.2%  cold 
water  fish  gelatin  (Sigma),  0.5%  BSA,  in  PBS).  Primary  antibodies 
recognizing  ATRX  and  DAXX  (Bethyl  Laboratories,  Inc.,  see 
above)  were  diluted  in  PBG  and  incubated  on  cells  for  1.5  h  at 
37°C.  Following  three  washes  with  PBG,  cells  were  incubated  with 
FITC-conjugated  donkey  anti-rabbit  secondary  antibody  (Life 
Technologies),  diluted  in  PBG,  for  30  min  at  37°C.  Cells  were 
washed  three  times  with  PBS  and  counterstained  with  4',6'- 
diamidino-2-phenylindole  (DAPI). 

Quantitative  RT-PCR 

RNA  was  extracted  from  cell  lines  using  the  RNeasy  kit 
(Qiagen).  0.5  pg  of  RNA  was  reverse  transcribed  using  RT- 
Advantage  kit  (Clontech)  in  a  volume  of  20  pi  according  to  the 
manufactures  instructions.  Following  reverse  transcription  180  pi 
of  water  was  added  and  5  pi  was  used  for  qPCR  with  specific 
primers  and  SYBR  green  mix  (Applied  Biosystems).  Primers: 
DAXX  5':  AGACGGTTTCTGAGCATCATC;  DAXX  3': 
AGAGGAGCTAGGGGCTTCTG;  TERT  5  b  GCCTTCAA- 
GAGCCACGTC;  TERT  3':  CCACGAACTGTCGCATGT. 

Analysis  of  copy  number  alterations 

Affymetrix  6.0  SNP  data  were  generated  at  the  Broad  Institute 
for  22  AET  cell  lines  and  6  single  cell  AET  (JECF-6/T.1R)  and 
non-AET  (HCT116)  cell  line  clones.  Single  cell  clones  were 
generated  by  fluorescence  activated  cell  sorting  (FACS)  sorting 
into  96  well  plates.  Normalized  copy  number  estimates  (log2 
ratios)  were  made  and  segmented  by  the  Circular  Binary 
Segmentation  algorithm  as  previously  described  [58].  The 
GISTIC  2.0  algorithm  was  performed  as  previously  described 
on  the  resulting  segmented  copy  number  data  from  the  22  ALT 
lines  [39,59].  The  boundaries  of  the  peak  amplified  and  deleted 
regions  identified  by  GISTIC  2.0  were  set  with  at  least  95% 
confidence  to  include  the  target  gene(s). 

Functional  analysis  of  DAXX  and  ATRX 

SV40-transformed  BJ  fibroblasts  at  PD  38  were  infected  with 
either  ATRX-  or  DTXX-specific  shRNAs  (Sigma)  with  the 
following  target  sequences:  ATRXsh589:  GCAGATTGATAT- 


GAGAGGAAT;  ATRXsh590:  CGACAGAAACTAACCCTG- 
TAA;  ATRXsh592:  CCGGTGGTGAACATTLAGAAAT;  DAX- 
Xsh800  :  GTLAGGGATGGACTAAGCTTLAT;  DAXXsh801: 
TCACCATCGTTACTGTCAGAA;  DAXXsh802:  GCCACA- 
CAATGCGATCCAGAA.  LacZ-specific  shRNA  was  used  as 
negative  control,  while  infection  with  hTERT  was  used  as  a 
positive  control.  Cells  were  selected  with  puromycin  (0.5  pg/ml) 
for  2  days,  and  at  least  1 0®  cells  of  each  sub-line  were  maintained 
in  culture  for  110  days.  Cultures  were  split  at  <80%  confluence, 
and  growth  medium  was  replaced  every  3  days.  Cells  were 
counted  at  each  passage  and  cumulative  PDs  were  recorded.  The 
suppression  of  ATRX  in  HeLal.3  cells  was  achieved  with 
ATRXsh592. 

T-SCE  assay  for  Telomeric  Recombination 

Cells  were  grown  in  10  pM  BrdU:BrdC  (3:1)  for  16  hr  with  the 
addition  of  0.1  pg/ml  colcemid  (Roche)  for  the  final  2  h.  Slides 
were  treated  with  0.5  mg/ml  RNAse  A  for  10  min  at  37°C, 
stained  with  0.5  pg/ml  Hoechst  33258  in  2x  SSC  for  15  min  at 
room  temperature,  and  exposed  to  365-nm  UV  light  (Stratalinker 
1800  UV  irradiator,  5400 J/m^).  Following  digestion  with 
Exo nuclease  III  (Promega,  10  U/pl,  2x10  min)  at  room  temper¬ 
ature,  slides  were  dehydrated  through  an  ethanol  series  (70%, 
90%,  100%)  and  incubated  sequentially  with  TAMRA-TelG  5'- 
[TTAGGG]3-3'  and  EITC-TelC  5'-[CCCT7kA]3-3'  probes  at 
room  temperature.  The  percentage  of  chromosome  ends  display¬ 
ing  telomeric  sister  chromatid  exchanges  (T-SCEs)  was  calculated 
from  four  independent  experiments. 

SKY 

Cells  were  treated  with  colcemid  for  1  hr,  harvested,  pelleted, 
re-suspended  in  a  hypotonic  solution  of  0.075  M  KCl  for  18  min, 
fixed  in  Carnoy’s  fixative  (3:1  methanohglacial  acetic  acid),  and 
washed  four  times  with  Carnoy’s  fixative.  All  fixed  samples  were 
spread  on  slides  for  staining  or  hybridization.  Chromosomes  were 
stained  with  Giemsa  or  Hoechst  to  visualize  chromosomal 
abnormalities.  Spectral  karyotyping  (SKY)  was  performed  on 
mitotic  samples  according  to  the  SkyPaint  DNA  kit  H-5  for  human 
chromosomes  procedure  (Applied  Spectral  Imaging,  SKY000029) 
and  imaged  on  a  Nikon  Eclipse  E6000  microscope  equipped  with 
the  SD300  Spectracube  and  Spectral  Imaging  acquisition 
software. 

G2/l\/l  checkpoint  analysis 

To  assess  G2/M  checkpoint  initiation,  subconfluent  cells  were 
exposed  to  y-irradiation  and  harvested  at  1  hr  with  non-irradiated 
cells  for  mitotic  index  measurement  by  flow  cytometry.  Briefly,  cell 
pellets  were  fixed  in  70%  ethanol,  washed  twice  with  cold  PBS, 
and  permeabilized  with  0.25%  Triton  X-100  in  PBS  for  15  min 
on  ice.  Cells  were  then  rinsed  with  PBS  containing  1  %  BSA  and 
stained  with  anti-phospho-histone  H3  SIO  antibody  (Cell  Signal¬ 
ing)  for  1.5  hr  at  room  temperature.  Following  two  washes  with 
PBS  containing  1%  BSA,  cells  were  stained  with  a  FITC- 
conjugated  donkey  anti-mouse  secondary  antibody  (Eife  Technol¬ 
ogies)  for  30  min  at  room  temperature.  Cells  were  washed  twice 
with  PBS  and  then  stained  with  propidium  iodide  (25  pg/ml)  in 
the  presence  of  0. 1  mg/ ml  RNase.  Plow  cytometry  was  performed 
on  a  BD  Biosciences  EACSCalibur,  and  the  percentage  of  mitotic 
cells  was  determined  as  those  that  were  EITC-positive  and 
contained  4N  DNA  content.  To  assess  maintenance  of  the  G2/M 
checkpoint,  nocodazole  (1  pg/ml)  was  added  1  hr  after  irradia¬ 
tion,  and  cells  were  ineubated  at  37°C  for  16  hr.  Cells  were  then 
harvested  and  processed  as  described  above.  To  account  for 
differences  in  the  rate  of  cell  cycle  progression  between  cell  lines, 
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the  mitotic  index  of  irradiated  cells  was  normalized  to  that  of  non- 
irradiated  cells  treated  with  nocodazole  for  16  hr. 

Micronucleation 

Asynchronously  growing  cells  were  fixed  and  examined  for  the 
presence  of  micronuclei  after  DAPI  staining  for  DNA.  Data  were 
obtained  from  three  independent  experiments  with  >  1 00  nuclei 
examined  in  each. 

C-circle  assay 

For  the  C-circle  assay,  DNA  from  cell  lines  was  isolated  using 
the  DNeasy  Blood  and  Tissue  Q,iagen  kit,  digested  with  Hinfl  and 
Rsal,  and  quantified  using  Hoechst  fluorimetric  analysis.  DNAs 
were  diluted  to  approximately  3  ng/pl,  quantified  again,  and  the 
dilutions  were  adjusted  to  exacdy  3  ng/pl.  The  assay  was 
performed  using  30  ng  of  quantified  DNA  as  described  by  Henson 
et  al.  [8]  using  dot-blotting  with  an  end-labeled  [AACCCT]4 
telomeric  oligonucleotide  probe.  C-circle  values  in  Table  1  were 
derived  from  three  independent  DNA  preparations  for  each  cell 
line. 

DSB  repair  assays 

Cells  grown  on  glass  coverslips  were  irradiated  with  0.5  Gy  IR 
and  incubated  at  37°C  for  1  hr  or  24  hr.  Irradiated  and  non- 
irradiated  cells  were  fixed  with  3%  paraformaldehyde  for  10  min 
at  room  temperature,  washed  three  times  with  PBS,  and 
permeabilized  with  0.5%  Triton  X-100  for  10  min  on  ice.  Cells 
were  then  washed  four  times  with  PBS  and  blocked  for  30  min  at 
room  temperature  with  PBG  (0.2%  cold  water  fish  gelatin  (Sigma), 
0.5%  BSA,  in  PBS).  Primary  53BP1  (Novus)  antibody  was  diluted 
in  PBG  and  incubated  on  cells  for  1  hr  at  37°C.  Following  three 
washes  with  PBG,  cells  were  incubated  with  FITC-conjugated 
donkey  anti-rabbit  secondary  antibodies  (Life  Technologies), 
diluted  in  PBG  for  30  min  at  37°C.  Cells  were  washed  three 
times  with  PBS  and  counterstained  with  4',6'-diamidino-2- 
phenylindole  (DAPI).  The  percentage  of  cells  containing  >10 
53BP1  foci  per  nucleus  was  calculated  from  three  independent 
experiments. 

For  the  PFGE  assay,  subconfluent  cells  were  irradiated  with  100 
Gy  IR  and  incubated  at  37°C  for  24  hr.  Irradiated  and  non- 
irradiated  cells  were  washed  once  with  PBS,  trypsinized,  and 
counted.  One  million  cells  per  sample  (in  duplicate)  were 
centrifuged  and  washed  a  second  time  with  PBS.  Cell  pellets 
were  resuspended  in  50  pi  10  mM  Tris  pH  7.2,  20  mM  NaCl, 
50  mM  EDTA  and  pre-warmed  at  55°C  for  10  min.  The  cell 
suspension  was  then  mixed  with  50  pi  of  2%  agarose/PBS  at  55°C 
and  solidified  in  plastic  molds  (BioRad).  Agarose  plugs  were 
incubated  with  1  mg/ml  Proteinase  K,  0.2%  sodium  deoxycho- 
late,  1%  N-lauroylsarcosine-sodium  salt,  100  mM  EDTA,  pH  8.0 
for  20  hr  at  50°C  and  then  washed  four  times  for  1  hr  with 
20  mM  Tris/HCl,  50  mM  EDTA,  pH  8.0.  Duplicate  plugs  were 
irradiated  with  600  Gy  IR,  and  all  samples  were  equilibrated  in 
0.5  X  TBE  for  30  min.  Agarose  plugs  were  inserted  into  the  wells 
of  a  1  %  agarose  gel  and  sealed  with  1  %  agarose.  PFGE  was 
carried  out  using  a  CHEF-II  apparatus  in  0.5  x  TBE  buffer  at 
6  V/cm  with  pulse  times  of  70  sec  (15  hr)  and  120  sec  (11  hr). 
Following  electrophoresis,  the  gel  was  stained  with  0.5  pg/ml 
ethidium  bromide  in  TBE  for  30  min. 

Supporting  Information 

Figure  SI  C-circle  assay  on  the  ALT  cell  line  panel.  Example  of 
dot-blot  results  of  C-circle  assay  on  ALT  lines  and  non-ALT 


negative  control  (BJ/hTERT/SV40,  blue).  Amounts  of  DNA  used 
for  the  assay  are  indicated. 

(TIF) 

Figure  S2  Lack  of  ALT-mediated  immortalization  upon  sup¬ 
pression  of  ATRX  or  DAXX.  A,  Q;RT  PCR  assay  for  the 
expression  of  hTERT  mRNA  in  the  indicated  cell  lines.  B,  C- 
circle  assay  on  the  indicated  cell  populations.  JFCF-6/T.  IF  is  a 
non-ALT  cell  line  that  serves  as  a  negative  control.  The  positive 
control  isJFCF-6/T.  1J/1-3C.  C,  Q;RT  PCR  assay  for  expression 
of  DAXX  mRNA  in  the  indicated  cell  populations.  D, 
Proliferation  of  SV40-transformed  BJ  fibroblasts  infected  with 
the  indicated  hTERT  or  shRNA  retroviruses. 

(TIF) 

Figure  S3  Immunoblotting  for  shelterin  components  in  ALT. 
Immunoblots  for  the  indicated  shelterin  components  in  whole  cell 
extracts  of  the  indicated  cell  lines. 

(TIF) 

Figure  S4  Analysis  of  TERRA  levels  in  ALT.  Northern  blots  of 
total  RNA  from  the  indicated  cell  lines  probed  for  TERRA  with  a 
C-strand  telomeric  repeat  probe.  The  ethidium  bromide  staining 
pattern  of  the  gel  is  shown  with  the  ribosomal  RNAs  indicated. 
The  bar  graphs  show  the  relative  expression  levels  of  TERRA 
derived  from  3-5  independent  experiments  and  standard 
deviations.  Telomerase-positive  controls,  BJ/hTERT/SV40  and 
HeLa,  are  shown  in  blue. 

(TIF) 

Figure  S5  Regions  of  recurrent  amplification  and  deletion  in  22 
ALT  lines.  GISTIC  2.0  analysis  of  deletion  (blue  lines,  left  panel) 
and  amplification  (red  lines,  right  panel)  events  identifies 
significandy  recurrent  peak  regions  (top  25  labeled  by  cytoband). 
False  discovery  rates  (q-values;  x-axis)  are  plotted  by  genomic 
position  (y-axis)  with  the  green  line  indicating  the  0.25  cut-off  for 
significance. 

(TIF) 

Figure  S6  DDR  signaling  in  ALT  lines.  Immunoblot  showing 
phosphorylation  of  CHK2  before  and  after  treatment  with  IR  in 
the  indicated  ALT  lines  and  the  RPE/hTERT  control  (blue). 
(TIF) 

Figure  S7  Example  of  assay  for  DSB  repair  kinetics.  The  graph 
shows  the  %  of  cells  with  >10  53BP1  foci  scored  by  IF  on  >100 
cells.  The  indicated  cell  lines  were  either  not  treated  with  IR, 
treated  with  0.5  Gy  and  incubated  for  1  hr  or  for  24  hr.  Data  from 
two  experiments  were  averaged.  Control  telomerase-positive  cell 
lines  (blue)  show  disappearance  of  IR-induced  foci.  Each  of  the 
ALT  lines  has  residual  IR-induced  DSBs  at  24  hr. 

(TIF) 

Table  SI  ATRX  exons  deleted  in  the  analyzed  ALT  lines. 
(XLSX) 

Table  S2  List  of  sequenced  genes. 

(XLS) 

Table  S3  Sequence  alterations  in  the  set  of  299  genes  analyzed 
in  the  ALT  cell  lines. 

(XLSX) 
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